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This  report  is  the  final  report  covering  work  performed  for  the 
Electro-Optical  Technology  Program  Office  at  the  Naval  Research  Laboratory 
under  Contract  N00173-77-C-0236.  The  purpose  of  this  contract  is  to  develop 
a system-level  computer  code  to  evaluate  the  performance  of  Infrared  Search 
and  Track  systems.  The  Infrared  Search  and  Track  System  Evaluation  Model 
(IRSTEM)  resulting  from  tnis  contract  is  a deterministic  code  which  models 
the  background  scene,  atmospheric  effects,  sensor  design  characteristics 
anu  al’  'ws  flexibility  of  incorporating  a variety  of  signal  processing 
options . 


This  report  summarizes  the  work  documented  in  the  three  quarterly 
reports  submitted  under  this  contract  and  provides  a summary  of  the  theory 

behind  the  models  used  In  IRSTEM.  a program  summary  and  a user's  manual 
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SECTION  1 
INTRODUCTION 

Wide  field-of-view  Infrared  Search  and  Track  (IRST)  systems  are  de- 
signed to  search  the  engagement  space  and  detect  incoming  threats  based  on 
the  IR  signatures  generated  by  the  exhaust  plume  or  aerodynamical ly  heated 
structure  (see  Figure  1).  These  systems  complement  RF  radar  because  they 
are  passive  and  because  they  operate  in  a different  portion  of  the  electro- 
magnetic spectrum  which  makes  it  difficult  to  jam  both  the  radar  and  the 
surveillance  sensor  simultaneously.  On  the  other  hand,  one  of  the  primary 
limitations  of  the  IR  surveillance  sensor  is  its  ability  to  detect  targets 
at  long  ranges  without  excessive  false  alarms  due  to  background  clutter. 
There  are  a variety  of  options  available  to  minimize  the  false  alarms  and 
increase  the  probability  of  detection  including  judicious  selection  of  the 
spectrum  bands,  focal  plane  arrays  and  signal  processing.  The  large 
number  of  options  make  it  impractical  to  experimentally  determine  the 
optimal  design,  but  there  has  been  a general  lack  of  analysis  tools 
available  to  complement  the  experimental  programs. 

The  purpose  of  the  IRST  evaluation  model  (IRSTEM  *s  to  provide 
the  analyst  with  a flexible  tool  for  determining  the  perron, ance  of  wide 
field-of-view  surveillance  sensors  under  realistic  battlefield  conditions. 
Flexibility  is  provided  by  allowing  the  user  to  specify  the  sensor  charac- 
teristics including  both  the  hardware  design  and  software  logic.  Realism 
is  provided  by  allowing  the  user  to  input  a spatially  and  spectrally 
varying  background  scene  which  can  include  clouds,  sun  glints  and  other 
potential  false  targets.  The  ability  of  the  search  set  to  to  then  detect 
the  target  can  be  analyzed  using  the  user's  specified  hardware  and  signal 
processing  since  the  code  predicts  the  time-dependent  signals  from  the 
individual  detector  preamplifiers. 

IRST  performance  analyses  generally  are  made  assuming  the  background 
noise  is  Gaussian  (in  magnitude)  and  white  (in  frequency).  These  assumptions 
greatly  simplify  the  computations  and  reduce  the  probability  of  detection 
(P^)  and  false  alarm  rate  (FAR)  computations  to  analytical  expressions: 
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Ficjure  1.  IRST  Schematic. 
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Hy  * irradiance  of  target  at  the  aperture  (w/ciir'/ 

T a threshold  level  for  clutter  rejection  (w/cm) 

* S/N  • NEI 

f \ , f-  * lower  and  upoer  video  cutoff  freguencies. 
respectivelv  (hr) 

on  * standard  deviation  of  background  or  system 
noise  (w/cnr) 

If  one  wishes  to  Incorporate  the  impact  of  spatial  frequency  variations  on 
system  performance  without  sacrificing  mathematical  simplicity,  the  power 
spectral  density  or  Wiener  spec tr-i.o , of  the  background  scene  can  be  used 
to  establish  the  magnitude  of  o . That  is. 


J P(K) 


where 


IFOV  of  the  detector  (sr) 
spatial  frequency 
1-D  Wiener  Spectrum 


a.b  » constants 

Using  such  a formulation  over  the  white  noise  approach  does  more  realisti- 
cally weight  the  average  power  in  a typical  background  scene  to  the  lower 
spatial  frequencies  but  it  Inadequately  models  the  high  frequency  spikes 
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which  tend  to  dominate  the  false  alarm  problems.  Hence,  the  Wiener  spectrum 
approach  tends  to  underestimate  the  false  alarm  rate. 

In  addition  to  the  basic  accuracy  problems,  the  statistical 
approach  to  IRST  analyses  also  suffer  from  the  problem  of  inflexibility.  The 
formulations  above  are  based  on  simple  threshold  crossing  logic.  And  since 
the  time-dependent  signal  response  from  each  pixel  is  lost  in  the  statis- 
tics, these  formulations  do  not  lend  themselves  to  evaluations  of  various 
signal  processing  logics  (e.g.  multi-color  discrimination,  spatial  discrimi- 
nation, etc.).  Therefore,  it  appeared  reasonable  to  embark  on  a new  Ped- 
ology which  is  embodied  in  IRSTEM. 

IRSTEM  is  structured  such  that  the  user  can  specify  each  background 
element  or  pixel  (typically  4096  elements),  each  of  which  can  differ  from 
its  neighboring  elements.  The  target  signature  and  potential  false  tar- 
gets which  can  vary  with  time  are  added  and  computation  of  the  signals 
arising  from  the  detectors/preamplifiers  is  made  based  on  the  spectrally 
varying  scene.  These  signals  are  modified  to  incorporate  the  spatially  or 
temporarily  varying  effects  (e.g.  the  optical  transfer  function  and  the 
electronic  response  function).  These  degradations  are  modelled  by  trans- 
forming the  signal  arrays  into  • » a.ency  space  and  then  modifying  the 
with  the  appropriate  transfer  functions  before  transforming  the  entire  arm- 
back  to  the  time/space  domain.  At  this  point,  the  signal  is  further  degraded 
by  random  system  noise  which  then  results  in  a time-dependent  signal  from 
each  detector 'preampl i fier  which  is  used  bv  the  signal  processing  routine. 

The  code  discussed  in  this  report  does  not  include  the  signal  processing 
options  (e.g.  spectral  and  spatial  discrimination,  signature  growth  and 
adaptive  threshold),  because  the  discrimination  techniques  and  the  order 
in  which  they  are  applied  depend  on  the  specific  system  design.  To  deter- 
mine the  range  at  which  the  sensor  can  discriminate  the  target  from  the 
background,  the  user  simply  manipulates  the  preamp  signals  in  exactlv  the 
same  manner  in  which  the  signals  are  processed  by  the  hardware  in  tne  field. 

This  report  summarizes  the  development  of  the  initial  IRSTEM  code 
and  includes  information  contained  in  the  three  quarterly  reports  submitted 
under  this  contract  (refs.  1,  2 and  3).  It  should  be  noted,  however,  that 
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some  of  the  information  included  in  those  reports  has  not  been  duplicated  here 
in  the  interest  of  brevity.  The  structure  of  the  report  includes  a theoreti- 
cal section  immediately  following  this  introduction.  In  that  section  two 
items  not  covered  in  previous  reports  (i.e.,  the  background  model  and  the 
sky/haze  model)  are  discussed  in  addition  to  the  logic  flow  and  oasic  method- 
ology for  each  of  the  major  code  elements.  Whenever  the  terms  are  not  de- 
fined for  the  equations  within  the  theory  section,  tiie  reader  should  refer  to 
the  List  of  Symbols  at  the  front  of  this  report.  Section  3 is  a concise  sum- 
mary of  the  codes  and  major  subroutines.  It  is  followed  by  a user's  manual 
including  input  and  output  instructions  and  a list  of  major  variables  used 
in  the  code  to  facilitate  code  operation  and  troubleshooting.  A sample  program 
is  provided  to  give  the  user  a checkpoint  to  ensure  the  proper  operation  of 
the  code,  and  ^ lir*'’'g  is  included  in  Apoendi  V 

It  sroi.  : >e  noted  tnat  some  of  the  features  me  joed  in  1RSTEM 
as  documented  iere  are  the  result  of  work  done  under  other  contracts. 

These  items  include  the  target  signature  and  the  target  signature  model 
which  was  done  under  SAI  Contract  N00039-77-C-0409  from  r.he  Naval  Elec- 
tronics Common.  Three  other  features  (the  variable  b ' ;round  array  speci- 
fication, inout/output  representation  in  wavelenctn  -ns*ead  of  wave  numbers, 
and  modification  of  the  integration  routine)  were  ,'unded  under  AF  Contract 
No.  F33615-78-C-151 1 . The  code  documented  in  this  report  is  only  the 
basic  program  which  is  expected  to  undergo  extensive  modification  and 
expansion  as  analyses  are  conducted  with  it.  A discussion  of  the  basic 
features  and  limitations  of  the  code  as  well  as  a list  of  recommendations 
for  potential  expansions  are  included  in  the  final  section  of  this  report. 


Section  2 

IRSTEM  LOGIC  AND  THEORETICAL  BASIS 


2.1  Logic  Flow 

The  basic  elements  of  IRSTEM  are  listed  in  Table  1 and  graphically 
depicted  in  Figure  2.  The  engagement  simulator  consists  of  the  user  speci- 
fied target  dynamics  and  the  specified  location  and  time  of  day  and  year 
for  the  sensor.  The  latter  specifications  are  used  to  determine  the  sun 
location  for  use  in  the  glint  and  solar  scattering  computations.  At 
present  the  target  dynamics  model  is  an  extremely  simple  one  which 
assumes  that  the  target  is  always  radially  in-bound  with  no  aspect  or  sig- 
nature change  during  the  engagement.  The  only  trajectory  inputs  required 
are  the  starting  position  of  the  target  (Rs),  its  velocity  (Vy)  and  the 
maximum  time  at  which  the  engagement  is  to  cease  (T  ).  However,  the 
engagement  simulator  is  structured  such  that  other  modes  of  operation 
(e.g.  crossing  targets,  targets  with  varying  flight  paths  and  targets 
with  changing  signatures)  can  be  incorporated  at  a later  time.  The 
geometries  associated  with  the  enqaqenent  space  are  discussed  in  sub- 
section 2.2  and  followed  by  discussions  of  the  background  models  and 
target  signature  models  in  the  following  subsections. 

At  present  the  false  target  information  is  assumed  to  be  caused 
by  sun  glints,  although  the  user  could  incorporate  flares  and  other 
decoys,  if  desired.  Sunglints  are  modelled  in  IRSTEM  as  randomly 
occurring  events  in  specific  cells  whose  location  is  specified  by  the 
user.  The  glint  irradiance  (or  voltage  response)  is  added  to  the 
aperture  irradiance  due  to  the  diffuse  solar  scatter  and  thermal  emission 
of  the  background.  Even  though  the  glint  areas  are  small  compared  to  an 
individual  background  pixel,  the  glint  irradiance  dominates  the  signal 
response  because  the  energy  propagates  as  a near-collimated  beam, whereas 
the  diffuse  solar  component  spreads  over  a hemispherical  solid  angle. 


Functional  Relationships  of  Primary  Code  Elements 


Table  1 

Basic  Elements  of  IRSTEM 


Engagement  Simulator 

• 

Background  (BKG)  Tape 

• 

Target  Signature  Tape  (TST) 

• 

False  Target  Tape  (FTT) 

• 

Propagation  Model  (L0WEM) 

• 

Sensor  Description  Input 

• 

Object  Fi ie 

0 

Image  File 

1 • 

Signal  Processing  Logic 

The  key  assumptions  in  the  sunglint  model  are: 

(1)  All  glint  areas,  Aq,  are  flat,  square  surfaces  50  cm  on  a 
side.  This  means  that  for  a background  range  of  10  km  the 
glint  area  is  about  6%  of  the  pixel  area  if  the  pixel  IFOV 
is  0.2  x 0.2  mr. 

(2)  The  solar  irradiance  incident  on  the  glint  surface  is  con- 
stant for  all  engagement  and  environmental  conditions.  The 
spectrum  used  in  IRSTEM  is  shown  in  Figure  3. 

(3)  The  reflected  solar  energy  propagates  to  the  aperture  with  a 
beamspread  which  matches  the  angle  subtended  by  the  sun 
(00=0.53  degrees) 

(4)  The  glint  surface's  reflectivity  is  0 or  0.8  depending  on 
the  luck  of  a series  of  random  draws. 

The  flint  irradiance,  H<.g(  ,R),  received  by  the  sensor  is 


hsgU>R) 
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The  propagation  effects  (i.e.  path  transmission  and  radiation)  are 
modelled  in  the  subroutine  called  LOWEM,  which  is  roughly  equivalent  to 
the  new  LOWTRAN  IV  atmosphtric  model.  It  also  includes  an  optional  water 
vapor  model  based  on  the  work  of  White  at  the  Atmospheric  Science  Labora- 
tory in  place  of  the  Burch  model.  To  utilize  this  routine,  the  standard 
LOWTRAN  inputs  (ref.  4)  are  required,  however,  one  of  the  planned 
improvements  is  to  reduce  the  propagation  subroutine  to  only  those 
LOWTRAN  features  which  are  required  for  IRSTEM.  This  reduction  should 
decrease  the  running  time  and  also  decrease  the  core  size  requirements. 

The  sensor  is  described  through  a series  of  inputs  for  the  spec- 
trally and  spatially  variant  parameters.  Spectrally  variant  parameters 
include  the  optical  transmission  function  (x  ) and  the  system  detectivity 
(D*).  The  detector  array  description,  optical  transmission  function, 
the  detector's  spatial  filter  function  and  the  electronics  filter 
function  are  examples  of  the  spatially  variant  effects.  The  latter 
three  transforms  are  discussed  in  more  detail  in  Section  2.6. 

The  key  simplifying  assumption  in  the  IRSTEM  methodology  is  that 
the  spectrally  dependent  and  spatial ly/temporal ly  dependent  effects  are 
separable,  '’his  assumption  is  graphically  depicted  in  Figure  4.  Note 
that  the  detector's  response  for  tne  spectral  band  of  interest 
(‘  3 1,  2,  3 * AAj,  AA>,  AAj  respectively)  is  computed  assuming  no 
image  degradation  from  atmospheric  turbulence,  platform  jitter,  optical 
train  distortion,  detector  shape,  or  processing  electronics.  These  effects 
are  incorporated  by  taking  their  modulation  transfer  function  (FL ) and 
multiplying  them  by  the  Fourier  transform  of  the  signal  voltage.  Note 
that  the  turbulence  and  jitter  effects  are  not  included  in  IRSTEM  at 
this  time. 

The  present  version  of  IRSTEM  stops  at  this  point  and  stores  the  time 
dependent  voltages  from  each  detector/preamp  in  the  system.  These  signals 
are  generated  for  each  passband  of  interest  as  the  detector  array  sweeps 
across  tie  background  scene  in  the  neighborhood  of  the  potential  target. 

Each  of  these  sweeps  are  updated  at  the  frame  time  of  the  detector,  which 
is  typically  of  the  order  of  1 to  2 seconds.  If  (a)  the  proper  background 
scene  has  been  inputted,  (b)  the  target  signature  is  realistic,  and  (c) 
the  atmospheric  inputs  and  sensor  inputs  included  in  the  transforms  are 
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ure  4 . Separation  of  Spatial  and  Spectral  Computations 


should  duplicate  those  voltages  which  one  could  measure  in  the  field. 

These  signals,  of  course,  are  not  sufficient  to  model  the  performance 
of  the  sensor  and  one  must  include  the  signal  processing  logic  in  addition 
to  the  hardware  and  phenomenological  descriptions.  For  example,  it  may 
oe  desirable  to  have  an  adaptive  threshold  which  requires  the  potential 
target  cell  to  be  several  times  the  average  signal  level  of  the  cells 
surrounding  it.  One  may  also  wish  to  require  that  the  signal  exceed  the 
threshold  level  for  several  scans  (e.g.  3 hits  out  of  5).  Or,  it  may 
be  desirable  to  compare  the  signal  ratios  between  spectral  bands.  One 
could  make  these  examinations  in  any  order  and  could  include  any  other- 
techniques  such  as  spatial  discrimination  or  target  signature  growth. 

Since  there  are  a wide  variety  of  options  available  to  the  user,  it 
was  decided  to  not  incorporate  a nominal  signal  processing  routine 
into  the  IRSTEM  code  at  this  time,  but  to  leave  the  signal  processing 


portion  of  the  program  to  the  user's  discretion.  Later  versions  of 
IRSTEM  may  include  a series  of  canned  programs  representing  specific 
hardware. 

2.Z  Geometry  Considerations 

The  geometric  computations  in  IRSTEM  are  somewhat  complicated 
because  three  different  reference  frames  are  used  in  the  program  (see 
Figure  5).  The  sensor  coordinate  system  (Figure  6)  is  oriented  such 
that  the  y-axis  always  coincides  with  the  centerline  of  the  background 
and  its  azimuth  argles  are  measured  counter-clockwise.  The  "true" 
reference  system  which  is  defined  by  the  zenith  angle  and  true  north, 
always  has  a fixed  orientation  in  space.  The  third  reference  frame 
is  defined  by  the  normal  to  the  cloud  surface  for  each  pixel  con- 
taining a cloud  element.  It  is  the  sun  and  observer  angle  with 
respect  to  this  normal  which  are  used  to  access  the  proper  scattering 
data.  The  orientation  of  these  reference  frames  and  the  pertinent 
angles  used  in  the  background  calculations  are  shown  in  Figure  7.  Note 
that  in  this  figure  RSP  refers  to  the  reciprocal  sun  position,  N is  the  cloud 
normal,  TN  is  true  north,  and  E is  East.  The  other  angles  are  defined  in 
the  list  of  symbols. 
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Figure  6.  Geometry  Sensor  Coordinate  System. 
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Since  the  background  scene  is  modelled  as  an  array  of  discrete 
pixels  high  by  J^x  across,  the  target's  angular  position  (0^,  (j>T) 
must  be  converted  to  scene  coordinates  (II,  JJ)  as  follows.  The  target 
position  is  compared  to  the  first  pixel  in  the  upper  left-hand  corner 
(Jlf  which  can  be  computed  knowing  the  centroid  of  the  scene 
(e^,  $ * 1.5707  rad)  and  the  pixel  size  (always  assumed  to  be  0.2  by 
0.2  mr).  That  is, 

rtr  - a l 

11  = 2 x'  Tor1*  + 1 


and 


= 5000(aT  - + IMAX *10~u ) + 1 


JJ 


4'*r  - 4*  \ 

2 x "10-V  * 1 


* 5000(tT  - 1.5707  ♦ JMftx*10“4)  ♦ 1 


Note  that  the  angular  coordinates  correspond  to  the  lower  right-*»and 
corner  of  the  pixel  (see  Figure  8).  Also  note  that  the  detector  array 
is  assumed  to  scan  the  scene  from  left  to  right  with  velocity,  v . 

2 . 3 Backgrounds  Model 

Five  thermal  backgrounds  have  been  generated  for  IRSTEM.  Each 
background  file  consists  of  coded  data  for  a 64  x 64  grid  of  0.2  x 0.2  mr 
pixels.  This  pixel  size  was  selected  to  represent  the  minimum  practical 
detector  IFOV. 

File  1 on  the  background  tape  consists  of  data  for  a partly  cloudy 
sky  with  a 23  km  haze.  The  field-of-view  consists  of  50%  sky,  and  50%  cloud 
at  a 10  km  range  (see  Figure  9).  The  cloud  is  made  up  of  a core  with 
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reflected  from  cloud. 
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Urban  scene  under  sky 
with  23  km  haze. 


Figure  9.  Five  Nominal  Backarounds. 


.5  km  thickness  and  -in  edge  with  .1  km  thickness;  Figure  9(a)  shows  the 
spatial  distribution  of  cloud,  cloud  edge,  and  sky  for  this  data  set. 
File  2 is  the  sea  viewed  toward  the  horizon  under  a night  sky.  The 
sea  is  assumed  tc  be  of  moderate  roughness,  and  the  sky  clear.  The 
spatial  distribution  of  sky  and  sea  is  as  shown  in  Figure  9(b). 

File  3 consists  of  the  sea  viewed  toward  the  horizon  with  a 
partly  cloudy  sky.  Here  the  sky  is  given  a 5 km  haze  and  radiation 
is  assumed  to  be  reflected  from  a cloud  of  10  km  thickness,  at  range 
10  km. 


File  4 consists  of  25%  land,  25%  sea,  and  50%  sky.  For  this 
background  a 23  km  haze  was  assumed,  with  radiation  reflected  from  a cloud 
of  10  km  thickness  at  a range  of  10  km. 


File  5 is  data  for  an  urban  scene  viewed  towards  the  horizon,  with 
sky  under  a 23  km  haze.  Data  used  was  for  Flint,  Michigan  and  represents 
a fairly  typicel  urban  distribution. 


Each  file  has  the  number  of  key  cards  in  the  first  line,  followed 
by  the  key  cards.  The  key  cards  are  formatted  (15,  4E10.3)  where  the 
integer  is  the  code  and  the  four  floating  point  numbers  give  the 
moaning  of  the  code,  in  the  following  order  temperature  (°K), 
emissivity  (%),  reflectivity  (%),  range  (meters).  This  information 
is  used  in  IRSTEM  to  generate  the  resulting  scene  radiance  and  corres- 
ponding signal  response.  These  backgrounds  are  not  intended  to  be  com- 
pletely representative  of  all  possible  spatial  variations  of  terrain  or 
sea  background.  However,  the  contrast  provided  by  mixing  scene  types 
in  a given  background  provides  a great  deal  of  the  variability  required 
for  sensor  performance  modelling. 


The  method  of  Fourier  synthesis  described  in  Rusbridae  (ref.  5) 
was  used  to  generate  background  scenes  with  proper  spatial  correlation 
between  scene  elements.  This  method  consists  basically  of  multiplying 
a two-dimensional  spatial  power  spectral  density  (PSD)  by  an  uncorv,''1ated 
random  array  and  taking  the  inverse  Fourier  transform.  The  resulting  array 
then  has  spatial  correlation  as  given  by  the  PSD.  One-dimensional  PDSs 
were  obtained  from  the  ERIM  (ref.  6)  and  Lockheed  (ref,  7)  field 
measurements.  An  example  temperature  contour  map  is  shown  in  Figure  10. 

In  general,  there  are  assumed  to  be  5 possible  temperatures , 2 possible 
ranges,  and  4 possible  emissivities  for  a total  of  40  different  -odes. 
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Figure  1CI  Moderate  Sea  Used  for  File  3 Ground  Temperature 
Mean  288.95°K  Std.  Dev.  .73°K 
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Selectivity  is  always  assumed  to  be  given  by  Kirchoff's  law 
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<2.4  Cloud/Haze  Model 

Clouds  appear  in  an  infinite  variety  of  sizes  and  shapes  and  are 
also  quite  variable  in  their  water  contents  and  distribtion  of  particle 
sizes.  During  the  daytime  a cloud  will  reflect,  transmit,  and  absorb 
incident  solar  and  sky  radiation  and  will  also  emit  thermal  radiation. 

A mathematical  model  can  be  developed  based  upon  the  radiative- transfer 
equation  which  w i < < predict  the  radiation  from  a cloud  provided  that  the 
microphysical  properties  are  known. 

The  major  obstacle  in  solving  the  radiati ve-transfer  equation 
involves  the  multiple  scattering  integral.  One  approximate  technique  is  that 
of  Turner's  (refs.  8,  . , 10)  model  which  uses  a modified  double-delta 
function  approximation.  The  phase  function  and  the  initial  solar  radiation 
fields  are  represented  as  delta  functions  in  angle  instead  of  the  highly 
forward  peaked  distribution  which  actually  exists.  Since  the  scattering  is 
very  pronounced  in  the  forward  direction  for  large  particles  characteristic 
of  clouds,  the  model  should  work  quite  well  for  cloud  radiation.  After 
considerable  calculation  involving  an  iterative  method,  a closed- 
formsolution  is  found  for  the  spectral  radiance  in  terms  of  the  scatter- 
ing phase  function. 

Rather  than  consider  all  possible  cloud  types,  IRSTEM  is  based  on 
the  tabulated  calculations  of  the  radiance  from  a cloud  element  for  a large 
variety  of  angles.  It  then  uses  that  element  as  a representative  facet  of 
the  cloud.  In  this  way  any  complicated  cloud  size  or  shape  can  be  simulated. 
We  must  then  consider  three  basic  regions  of  the  cloud.  Region  1 represents 
sunlight  reflected  from  an  ontlcally  thick  cloud;  region  2 represents  radia- 
tion which  is  transmitted  through  an  optically  thin  portion  of  the  cloud 
edge;  and  region  3 represents  radiation  which  is  emitted  from  the  optically 
thick  cloud.  In  addition,  the  radiance  for  a 5 Km.  and  a 23  Km.  visual 
range  hazes  has  been  calculated  and  tabulated  for  IRSTEM  by  Turner.  It 
was  a maritime  haze  with  a real  particulate  refractive  index  of  1.5. 


To  incorporate  fhe  cloud/haze  radiance  data,  the  user  first  iden- 


tifies the  background  cells  in  which  the  cloud  is  located  and  specifies 
the  direction  of  tne  norma i to  the  cloud  surface  for  that  cell.  The  code 
then  takes  the  appropriate  radiance  data  from  the  cloud/haze  data  files. 
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interpolating  points  if  required,  based  on  the  relationship  between  the  solar 
angles,  viewing  angle,  and  direction  of  the  surface  normal.  Graphical  repre- 
sentations of  portions  of  the  data  base  incorporated  in  IRSTEM  are  snown  in 
Figures  11  and  12  for  \=2  um.  All  units  are  given  in  mw/sr/cnr/um. 

2 . 5 Target  Signature  Model 

A oetailed  target  signature  model  has  been  incorporated  iniu  tne 
program  based  on  recent  SAI  work  under  NAVELEX  Contract  No.  N00039-77-C-040'J>. 
The  various  signature  components  are  treated  separately  to  give  a physically 
more  realistic  model.  The  radiant  intensity  from  the  rocket's  plume,  tail- 
pipe, and  aerodynamical ly  heated  skin  are  combined  to  give  the  overall  spec- 
tral radiant  intensity  as  a function  of  aspect  angle  9,  as  shown  in  equation  5. 
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Spectral  radiant  intensity  (w/sr-um) 

Plume  transmission 

Blackbody  spectral  radiance  (w/cm--sr-um) 
Projected  area  viewed  from  aspect  angle  Z (cm) 
Aspect  dependence  of  the  plume  emission 
Temperature  (°K) 
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p - plume 
E - exhaust  nozzle 
N - nose  of  missile 
S - skin  of  missile 

The  geometric  representation  of  a generic  missile  is  .shown  in  Figure  13. 

The  first  term  in  equation  5 is  the  plume  radiance,  which  generally 
dominates  the  signature.  The  second  term  is  the  tailpipe's  blackbody  radia- 
tion which  is  observable  only  at  large  aspect  angles,  and  the  last  two  terms 
odei  tr.e  skin  emissions.  It  should  be  noted  that  the  resulting  signature 
is  a function  of  \ and  0 only.  The  change  in  radiant  intensity  over  the 
projected  area  of  the  target  is  not  of  interest  since  we  are  interested  in 
analyzing  the  initial  detection  capability  of  the  IRST,  and  at  these  ranges 
the  target  appears  as  a point  source. 

The  plume  intensity  Jp  has  been  computed  using  the  ASDIR  signature 
code  for  a missile  at  a 20°  aspect  angle  along  with  a function  fp(e)  to 
;'•?  these  results  to  other  viewing  angles.  This  assumed  aspect  depen- 
dence is  shown  in  Figure  14. 

The  emission  from  the  heated  tailpipe  of  the  threat  is  represented 
by  the  second  term  in  equation  5.  That  term  is  the  blackbody  radiance 
function  inultipled  by  an  effective  area  for  the  tailpipe  and  a plume 
transmission.  The  effective  projected  area  varies  as  a function  of  viewing 
aspect  and  the  missile  type.  A zero  projected  nozzle  area  is  assumed  for 
forward  aspect  angles  (0  < 90°). 

For  purposes  of  calculating  the  tailpipe  emissions,  we  have  assumed 
that  the  plume  transmission  is  independent  of  aspect  angle.  The  plume  trans 
mission  function  used  is  based  (arbitrarily)  on  a 160°  aspect  angle.  This 
function  is  shown  in  Figure  15  along  with  a 110°  curve  for  comparison  which 
shows  that  the  assumption  to  not  vary  the  plume  transmission  as  a function 
of  aspect  angle  is  a reasonable  one. 
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Figure  15-  Transmission  through  the  Plume  for  a 110°  and  160 
Viewing  Geometry  to  the  Exhaust  Nozzzle. 


The  missile  nose  and  skin  emissions  are  computed  based  simply  on 
♦■he  oroduct  of  the  blackbody  function  and  a projected  area.  The  radiance 
from  aerodynamic  heating  is  divided  into  two  terms  so  that  higher  tempera- 
ture 'at  or  near  the  stagnation  temperature)  can  be  assumed  for  the  nose  of 
■ e • ,i'e  and  more  moderate  temperature  assumed  for  the  rest  of  the  skin. 

The  various  portions  of  a typical  missile  signature  and  the  result- 
ant total  signature  are  shown  in  Figures  16,  17,  18  and  19.  Note  that  these 
example  plots  are  based  on  viewing  the  missile  from  the  near,  which  is  not 
consistent  with  the  present  trajectory  mode  of  only  radially  in-bound  tar- 
gets (ref.  Section  2.1).  This  generalized  signature  model  which  includes 
the  tailpipe  emissions  has  been  developed  in  anticipation  of  a more  general  tra- 
jectory model  in  later  versions  of  I RSTEM  which  will  model  : rousing  targets. 

2 . c Sensor  Model 
2.6.1  Methodology. 

If  we  assume  the  projection  of  the  detector  array  moves  across 
the  background  scene  as  shown  in  Figure  20,  the  ideal  voltage  generated 
at  the  output  of  the  detector  due  to  any  given  pixel 
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assuming  0.2  x 0.2  mr  pixels. 

Equation  (6)  is  a manifestation  of  one  of  the  key  simplifying 
assumptions  used  in  IRSTEM;  namely,  the  separability  of  the  spectrally  and 
spatially  dependent  effects. 

To  compute  the  spatial ly-dependent  effects  (or  temporally  dependent 
effects  since  x = v$t),  we  consider  3 distortions:  the  sensor  optics,  the 
finite  detector  size,  and  the  electronics  response.  The  optical  distortions 
are  incorporated  by  convolving  the  ideal  voltage  with  the  optical  point 
!Ad  function  (PSF)  to  produce  an  effective  signal  response,  V 
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Figure  19-  Combined  Signature  for  a Typical  Missile 
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where 


o * ix  + jy 

To  convert  the  static  background  scene  into  the  time  varying  signals  arising 
from  the  n individual  detectors,  we  integrate  (or  sum)  the  background 
signals  arising  from  each  pixel  as  weighted  by  the  detector's  aoerture 
function. 


where 


(x,y)d(x  - v$t,y)  dxdy 


d(x,y)  = 


Lx  Lv 

1,  for  J x | <_  and  Jy  I <_ 

0,  elsewhere 


(8) 


Finally,  we  obtain  the  time  varying  signals  from  each  detector/ 
preamplifier  circuit  by  convolving  the  detector's  output  signal  with  the 
electronic  response  function 


vn(t) 


= v’(t)*hE(t) 

= j v;(f)hE(t  - t;)dt' 


(9) 


The  computations  could  be  made  directly  as  indicated  by  the 
equations;  however,  the  appearance  of  convolutions  immediately  suggests 
the  use  of  Fourier  transforms.  This  approach  is,  in  fact,  quite  compatible 
with  a digital  implementation  as  many  fast  discrete  transform  algorithms 
like  the  Fast  Fourier  Transform  or  FFT  are  now  available. 
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Defining  the  one-  and  two-dimensional  transform  pairs  as 


f(t)  = ^ J F(^)elL 
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we  can  rewrite  equations  7,  8 and  9. 
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vs  f Veff(vs  ’ ky)Hdet(vs  ,ky)dky 


and 


Vn( u)  = He(w)V^(m) 


Therefore  we  can  combine  equations  10,  11,  and  12  and  compute  v ( t ) as 


Vn(t)  = h ffh deal<kx-ky)Hopt(kx-ky)Hdet<kx'ky>Helec(Vx> 


ik  v t 


dk  dk 
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In  other  words,  the  transformed  ideal  voltage  array  is  modified  by  the 
appropriate  transform  for  the  optics,  detectors  and  electronics  to  create 
a new  array  which  is  inverse  transformed  to  generate  the  resultant  time- 
dependent  signal  response  from  each  detector. 


2.f.2  Sensor  Transforms. 

The  transfer  functions  which  are  currently  built  into  IRSTEM  are 
summarized  as  follows: 


Optics  - Circular  Clear  Aperture 
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(12) 


is  the  linear  extent 
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where 

HDet  = Detector  response 

A . = Detector  area  = L L 

u ^ y 

k = Spatial  frequency  vector,  lk  + jk 

x y 

sinc(x)  * sin(x)/x 

Electronics  - Low-Pass  Filter 


HElec(kxvs) 


whe  re 


v$  = Image  -can  rate  across  the  detectors 

fo  = Filter  cutoff  frequency 

F 3 Focal  length 

Ts  = Scan  time 


Even  though  these  transforms  are  currently  hardwired  into  IRSTEM. 
the  user  can  replace  them  with  alternative  expressions.  With  some  slight 
code  modification,  the  user  can  even  specify  the  physical  characteristics 
of  the  hardware  (e.g.  aperture  shape,  obscurations,  and  distortions)  and 
then  use  the  Fast-Fourier  Transform  to  compute  the  frequency  response 
function  when  an  analytics7  axpression  is  not  readily  available. 
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2.6.3  Other  Sensor  Considerations 

Selecting  Appropriate  Detector  Response 

The  image  array  which  results  from  the  inverse  transform  operation 
discussed  above  is  based  on  the  full  number  of  pixels  in  the  background 
scene.  Since  the  instantaneous  field-of-view  (IFOV)  of  the  detectors 
are  generally  larger  than  0.2  x 0.2  mr  pixels,  the  rows  corresponding 
to  the  center  of  the  detector  are  selected  to  represent  the  output 
signals  and  the  rest  of  the  array  is  discarded.  The  determination  of 
which  rows  correspond  to  detector  centers  proceed  as  shown  in  Figure  21. 
The  number  of  detectors  falling  fully  within  the  vertical  scene  dimension 
is  computed  using  the  detector  length  and  degree  of  overlap  with  the  next 
detector. 


The  detector  centers  are  at 

yi  a 2^  + (t-'lKLy-OVL) 

or  pixel  numbers 

II  = y./(2  x 10"*F)  + 1 

since  the  top  edge  of  the  first  detector  always  coincides  with  the  top 

edge  of  the  scene.  In  the  example  shown  in  Figure  21,  all  dimensions 
are  normalized  to  the  pixel  size  as  follows. 

Y - 18.5  OVL  = 1 

Ly  ’ 5 

which  implies  N * 4 with  detector  centers  at  3,  7,  11  and  15. 

System  Noise. 

Random  system  noise  (e.g.  thermal,  shot,  or  1/f  noise)  is  modelled 
in  IRSTEM  by  modifying  the  image  file  voltages  by  a random  number,  vn. 


Figure  21. 


Deter^ina-- ;n  Detector  Centers. 


» 


» 


h’ 


wh-'ch  is  uncorrelated  from  cell  to  cell  (see  Figure  22).  A random  number 
sequence  of  299  numbers  with  zero  mean  and  unity  standard  deviation  is  built 
into  the  code.  The  user  must  input  the  rms  noise  level  (or  standard  devia- 
tion since  we  are  assuming  Gaussian  statistics)  to  allow  the  program  to 
scale  the  random  numbers  properly. 

A simple  way  to  estimate  this  noise  level  is  to  compute  the  rms 
noise,  vn  as  follows: 
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where 


HE(f) 
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D* 


Load  resistance 

Electronic  filter  response  function 
Detector  noise  voltaae 

R j 

Responsi vi ty 
Detectivity 
Detector  area 


Note  that  the  system  noise  is  simply  the  detector  noise  modified  by  the 
filter  response  function.  This  response  function  must  be  integrated  over 
the  frequency  interval  of  interest  to  properly  account  for  the  additional 
noise  input  due  to  the  filter's  limited  bandpass  since  the  filter  response 
function  already  has  been  convolved  with  the  object  file,  the  optical 
point-spread  function,  and  the  detector  response.  If  we  have  an  ideal  low- 
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"'gure  22.  Random  Noise  Input  to  IRSTEM 


pass  filter  with  cutoff  f0  , the  integral  in  equation  (13)  is  simply  f 0 ; 
if  the  response  function  is  l/(l+f/f0),  the  result  is  irfn/Z.  This 
implies  the  rms  noise  input  can  be  given  by 


/V 


W peak  ’ V» 


where  K = 1 for  an  ideal  low-pass  filter 

= v'7/2  for  a more  realistic  low-pass  filter 


The  peak  values  for  R^  and  D*  are  shown  because  they  are  readily  available 
from  the  detector  manufacturers.  These  values  are  also  required  inputs  along 
with  D*(\)  to  compute  the  spectral  responsi vi ty.  That  is,  the  code  computes 
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peak 


We  do  not  require  Rd(A)  to  be  user  specified  since  this  information  is  not  as 


readily  available  as  R.  , D*  , , and  0*0) 

a peak  peak 


2.  7 Surra;-  Equat~,  ,:is 

The  basic  equations  used  in  IRSTEM  are  summarized  below  in  approxi- 
mate order  of  apcearance.  The  list  of  symDOis  at  the  front  of  this  report 
defines  each  parameter. 
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In  this  section  the  IRSTEM  code  is  reviewed  to  identify  its  major 
elements  and  the  methodology  for  combining  these  elements.  The  sub- 
routines are  listed  along  with  a functional  description  and  the  inputs 
and  the  outputs  associated  with  each  subroutine.  Specific  instructions 
for  operating  the  code  are  discussed  in  Section  4. 


3.1  Overall  Program  Flow. 

The  top  level  flow  diagram  for  IRSTEM  is  shown  in  Figure  23-  The 
program  begins  by  reading  tne  background  codes  associated  with  each  pixel, 
'hen  the  code  keys  are  read  to  specify  the  temperature  emissivity,  reflec- 
tivity and  range  associated  with  each  unique  background  code.  Since  these 
codes  may  not  be  numbered  consecutively,  the  codes  are  reordered  to  make  them 
sequential  reserving  codes  1 through  6 for  sky  and  cloud  data. 

After  the  background  information  has  been  fixed,  the  scene  is 
scanned  to  identify  duplicate  ranges  to  minimize  calls  to  LOWEM  later  in 
the  program.  Also  at  this  time,  the  minimum  and  the  maximum  ranges  for  the 
entire  background  scene  are  identified  although  they  are  not  used  in  this 
version  of  IRSTEM.  Anticipated  modifications  to  minimize  the  computational 
time  probably  will  use  this  minimum/maximum  range. 


Since  LOWTRAN  (i.e.  Subroutine  LOWEM)  requires  regular  wavenumber 
spacings  (Av  or  CV)  that  are  multiples  of  5 cm"1,  the  integrations  to 
form  tne  arrays  must  be  done  in  terms  of  \ spacings  which  correspond 

to  equal  v intervals.  These  spacings  are  fixed  in  the  early  part  of  the 
program  as  follows: 
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Figure  23.  IRSTEM  Flow  Chart. 


where  { } implies  integer.  The  integration  is  then  done  at  uneven  \ 
spacings  where  the  AA  interval  at  the  lower  wavelength  bound  is  equal 
to  the  spacing  specified  by  the  user. 

An  example  spacing  is  shown  in  Figure  24  for  the  4.42  - 4.5  pm 
passband  with  AA  =0.02  um. 

Av  = j 2. 05 1 5 = 10  cm* ' 

vL  = | 444 | 5 = 2220  cm'1 

u = 1 453.49|5  = 2265  cm'1 

The  upper  line  shows  t.ie  interval  of  interest  .vnic::  is  approximated  bv  the 
smallest  wavenumber  band  divisible  by  Av  (i.e.,  \,  - v ).  The  intervals 
for  Av  = 10  cm  spacing  begin  at  vL  = 2220  and  go  to  2270  cm  to  encompass 
v = 2265  cm  ‘.  The  corresponding  wavelength  interval  is  4.4053  - 4.505  um. 
This  particular  example  is  a worst  case  approximation  since  the  interval 
spacing  is  large  compared  to  the  total  passband  but  it  illustrates  metho- 
dology used  in  IRSTEM. 

Following  specification  of  the  intervals  to  be  used  in  the  integra- 
tion routines,  the  plume  intensity,  plume  transmission,  optical  transmission 
system  detectivity  and  solar  radiance  are  read  as  a function  of  wavelength 
and  interpolated  to  match  their  values  to  the  values  required  for  the  inte- 
gration routines.  These  operations  are  handled  in  the  subroutine  called 
INTER.  The  final  major  set  of  data  input  is  the  scenario  information 
including  both  the  target  trajectory  data  and  the  sun  location  information 
used  in  the  glint  calculations.  This  input  phase  of  the  program  is  then 
followed  by  a series  of  preliminary  calculations  before  entering  the  first 
major  due  loop  of  the  program. 

The  MXC0DE  loop  shown  in  Figure  23  creates  the  initial  background 
voltages  for  each  pixel  in  the  background.  This  computation  is  done  for 
each  of  the  static  background  inputs,  which  were  read  at  the  beginning 


Figure  24. 
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of  the  program,  but  it  bypasses  the  first  6 codes  which  are  associated  with 
the  sky  and  the  cloud  conditions.  The  loop  operates  as  follows:  A check 
is  made  to  see  if  the  range  to  the  background  cell  has  been  used  in  pre- 
vious calculations.  If  the  range  has  not  been  encountered  previously, 
the  voltage  due  to  the  path  radiance  and  the  diffuse  solar  scatter  are 
combined  to  form  V12.  These  voltages  represent  the  signal  levels  at  the 
output  of  the  detectors  in  the  absence  of  optical  distortions,  detector 
spatial  filtering  and  the  electronics  filtering.  These  calculations  are 
made  by  integrating  over  the  wavelength  bands  of  interest  (up  to  3 passbands) 
If  the  range  for  the  particular  background  code  has  been  used  previously,  the 
V i 2 and  path  transmission  calculations  are  bypassed  and  the  program  proceeds 
directly  to  the  computation  of  the  third  background  component  (the  voltage 
array  due  to  the  background's  thermal  emission).  Finally,  all  three  com- 
ponents are  combined  to  form  the  voltage,  which  uniquely  corresponds 

to  each  of  the  background  codes.  Once  the  voltages  for  all  of  the  codes 
have  been  computed,  the  voltages  are  assigned  to  the  appropriate  background 

cel1s  (Veff ^ ’ 

The  sky/cloud  subroutine  is  called  next  to  compute  the  signal  responses 
arising  from  either  clear  sky,  hazy  sky  or  clouds.  This  information  is  added 
to  the  Vet-^  array  and  the  entire  array  is  stored  in  case  the  user  wishes 
to  bypass  the  calculations  up  to  this  point  in  future  analyses.  The  cal- 
culations to  this  point  need  be  done  only  once  for  a given  run  of  IRSTEM. 

At  this  point  program  control  is  transferred  to  the  second  major  loop  - 
the  time  loop. 

The  time  loop  begins  by  calling  subroutine  GLINT,  which  modifies 
the  background  voltage  array  by  adding  signal  responses  due  to  glints 
to  the  cells  which  contain  potential  glint  points.  Whether  a glint  occurs 
on  a given  scan  is  dependent  on  the  luck  of  a random  draw.  The  trajectory 
routine  is  then  called  to  update  the  position  of  the  target  and  the  time 
clock.  Next,  the  target  signature  subroutine  is  called  to  compute  the 
signature  of  the  target  at  this  point  in  time  and  to  modify  the  voltage 
array  for  the  cell  containing  the  target.  If  the  target  exceeds  the  pixel 
size,  the  program  stops  since  the  code  is  currently  designed  to  handle 
only  point-source  targets. 


48 


Now  the  scene  array  is  ready  for  transformation  to  the  frequency 
domain  so  that  the  array  can  be  modified  by  the  sensor  induced  distortions. 
This  transformation  is  accomplished  by  railing  subroutine  FIX  which  con- 
trols the  calls  to  the  fast-Fourier  *.*ansform  subroutine  (FOURT).  Once 
the  transformation  is  complete,  the  transformed  array  is  modified  by  the 
optical  transmission  function,  the  detector  transmission  function  and  the 
electronics  transmission  function  to  create  a new  array  which  is  transformed 
back  to  the  time/space  domain.  Next,  random  noise  is  then  added  to  each  of 
the  scene  cells  based  on  the  normalized  rms  noise  input  specified  by  the  user 
and  the  result  of  a random  number  draw.  Finally,  those  rows  corresponding 
to  the  sensor  specified  by  the  detector  geometry  are  selected  for  manipula- 
tion by  the  signal  processing  routine.  The  rest  of  the  information  is 
dropped  from  the  file. 


At  this  point,  the  signal  processing  subroutine  could  be  called  by 
tne  user  to  analyze  the  performance  of  the  search  set.  As  noted  earlier, 
this  routine  is  not  included  at  this  time  because  of  its  unique  relationship 
to  each  specific  hardware  design.  The  data  presently  is  stored  on  a file  in 
anticipation  that  a signal  processing  routine  will  be  called  after  the  time 
loop  is  completed;  perhaps  the  : jnal  processor  will  even  be  written  as  a 
separate  program.  It  is  anticipated  that  the  signal  processing  routine  will 
analyze  the  raw  video  signals  generated  by  IRSTEM  to  determine  at  what  point 
the  target  can  be  detected  based  on  the  discrimination  techniques  employed. 


Subroutines  and  Functions 


The  subroutines  and  functions  called  by  the  main  program  are  listed 
alphabetically  in  Table  2.  Most  of  these  routines  are  discussed  further 
in  the  paragraphs  below.  Some  of  the  more  simple  routines  (e.g.,  Dl'MR  and 
INTEG)  are  not  reviewed;  the  LOWEM  subroutines  (HN03,  POINT  and  ANGL)  are 
also  not  reviewed  since  they  are  described  in  reference  4.  Summary  tables, 
for  each  subroutine  except  those  associated  with  LOWEM,  are  provided  at  the 
end  of  this  subsection  to  delineate  input, output  variables. 


Table  2. 

IRSTEM  Subroutines  and  Functions 


Subroutine 

ANGL 

DUMR 

FIND 

FIX 

FOURT 

GETANG 

GLINT 

HN03 

INTEG 

INTER 

LOWEM  . 

POINT 

SCALR 

SKYCLD 


Purpose 

LOWEM  subroutine  to  compute  zenith  angle. 

Dummy  read  routine. 

Selects  data  sets  for  sky/cloud  radiance 
interpolations. 

Calls  Fast  Fourier  Transform  and  reorders 
the  result  of  the  transformation. 

Fast  Fourier  Transform. 

Calculates  sun  and  observer  angles  with 
respect  to  the  cloud  normal. 

Calculates  signal  response  due  to  solar 
glints. 

LOWEM  Subroutine  for  HN03  extinction 
coefficients. 

Integration  routine  (trapezoidal  rule). 

Reads  Jp,  xp,  xQ,  and  H.  as  a function 

of  \ and  interpolates  values  needed  for 
integration.  Also  interpolates  cloud  radiance. 

Path  transmission  and  radiance  model 
(LOWTRAN  IV). 


LOWEM  subroutine  to  interpolate  absorber 
concentrations  as  a function  of  altitude. 

Scales  a^ray  for  printing. 

Calculates  voltage  response  for  sky  and 
cloud  background  elements. 


Table  2. 


(Continued) 


Subroutine 

Purpose 

SUNLOC 

Calculates  sun  location. 

TGTSIG 

Calculates  target  signature. 

TRAJR 

Calculates  target  position  as  a function 
of  time. 

Function 

BLKBDY 

Computes  the  blackbody  emittance  for  a 
specified  wavelength  and  temperature. 

INTERP 

Interpolates  plume  signature  angular 
dependence. 

NODD 

Generates  +1  or  -1  depending  on  whether 
the  function  is  even  or  odd. 

ii 
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Subroutine  FIND.  The  sky/cloud  spectral  radiance  data  is  stored  in 
5 files  (see  discussion  of  Subroutine  SKYCLD  and  the  model  methodology  in 
Section  2.4)  as  a function  of  the  solar  and  view  angles.  The  angles  for  the 
reflected  solar  and  haze  situations  are: 


solar  angle  (e0) 
relative  azimuth  angle  (<j>) 
view  angle  (e) , or 
supplement  view  angle  (o') 


0,  15,  30,  45,  60,  75,  85 
0,  45,  90,  135,  180 

0,  15,  30,  45,  60,  75,  85 


If  we  designate  each  spectral  radiance  data  file  (,\  = 2,  3,  4,  5,  6,  8 
10  and  14  urn)  with  a 3 parameter  code  corresponding  to  ( 9 0 , <J>,  0),  the  data 
set  (1,  3,  2)  implies  90  = 0°,  <j>  = 90°  and  0 = 15°.  If  90  = 15°  and  $ = 150° 
and  9 = 50°,  the  desired  data  would  lie  somewhere  between  the  (1,  4,  4)  and 
(1,  5,  5)  curves.  This  subroutine  calculates  the  ratios  between  the  two 
bounding  curves  through  three  successive  calls  from  SKYCLD. 

The  angles  for  files  2 and  3 (transmission  through  0.1  and  0.5  km 
thick  clouds  respectively)  are: 

azimuth  angle  (4>)  : 0,  45,  90,  135,  180 

view  angle  (01)  : 0,  2,  4,  6,  8,  10 

Here  the  solar  angle  (90)  is  not  explicitly  used  since  the  1 ine-of-sight 
between  the  observer  and  the  sun  is  the  reference  angle  (see  discussion  of 
subroutine  SKYCLD). 

Subroutine  FIX.  This  routine  calls  the  fast  Fourier  transform  and 
then  reorders  the  transformed  array.  This  re-ordering  is  necessary  as  the 
transform  has  been  formulated  for  values  on  the  interval  (0,  N-l)  which 
yield  2 half  periods  back  to  back.  It  is  necessary  to  move  the  origin  to 
N/2  as  shown  in  Figure  25  to  ensure  the  peak  of  the  transformed  spectrum 
corresponds  with  the  d.c.  component  of  the  transfer  functions. 
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Figure  25.  Shifting  the  Transformed  Signal. 


Subroutine  FOURT.  This  subroutine  performs  the  fast  Fourier  trans- 
form. The  subroutine  is  being  u-ed  with  no  modification.  FOURT  performs  an 
n-dimensiona 1 fast  Fourier  transform  on  an  n-dimens ional  array  of  complex 
data.  The  transform  performed  may  be  expressed  as  follows:  TRANSFORM 
(Jl,  J2,  ...)  = SUM  (DATA  (II,  12,  . . . )*W1**( ( II -1 )*( Jl-1 ) )*W2**( ( 12-1 )* 

( J2- 1 ))*...)  where  II  and  Jl  run  fr~m  I to  the  maximum  of  the  first 
dimension,  12  and  J2  run  from  1 to  the  maximum  of  tne  second  dimension. 

The  basis  of  the  routine  was  develooed  by  Cooly  and  Tukey.  The  FORTRAN 
code  being  used  was  developed  in  1965  at  Sandia  Corporation.  The  input 
data  for  the  transform  is  assumed  to  be  in  the  time  domain  and  equally 
spaced  at  intervals  of  OT.  The  transformed  values  are  considered  to  be  in 
frequency  space.  The  frequency  intervals  are  equally  spaced  from  0 to 
2 (N(I)-1)/(N(I)+DT)  where  N(I)  is  the  dimension  of  the  data  (64)  in  each 
of  the  I dimensions  (2).  FOURT  does  not  require  the  array  to  be  a multiple 
of  2,  but  the  computations  are  faster  if  this  criteria  is  met.  This  rou- 
tine was  chosen  because  of  its  flexibility  and  computational  speed.  A 
two-dimensional  transform  of  a 64  x 64  array  requires  less  than  one  second 
on  the  CDC  6600. 

Subroutine  GETANG.  Once  the  sun's  location  has  been  established, 
subroutine  GETANG  then  converts  this  information  along  with  the  observer’s 
position  into  the  cloud  reference  frame  before  calling  the  solar  scatter 
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data  table.  The  required  angles  are  30,  o and  o which  are  computed  using  the 
law  of  cosines  for  the  angles  projected  on  a celestial  sphere  (see  Ref.  3 
and  Section  2.7) . 

Subroutine  GLINT.  Immediately  upon  entering  the  time  loop  in  tne 
main  program,  the  glint  subroutine  is  called  to  modify  the  background 
voltage  array  for  the  cells  containing  potential  glint  points.  The  flow 
diagram  for  this  routine  is  shown  in  Figure  26.  If  this  is  the  first  time  the 
subroutine  has  been  called,  it  reads  in  the  number  of  glint  locations,  the 
location  of  potential  glints,  and  the  probability  of  glint  occurrence  and 
then  computes  the  corresponding  signals  at  the  output  of  the  pre-amps  due 
to  each  glint.  If  this  is  not  the  first  time  that  the  subroutine  has  been 
called,  the  signal  levels  computations  are  not  repeated,  and  program 
control  is  transferred  directly  to  the  loop  which  calculates  the  random  number 
for  each  glint  point.  If  this  random  number  exceeds  the  probability 
of  glint  occurrence  the  voltage  for  the  cell  under  examination  is  modified 
by  the  precalculated  glint  voltage  before  program  control  is  transferred 
back  to  the  main  program. 

Subroutine  INTER.  This  subroutine  serves  a dual  purpose  of  reading 
the  wavelength  dependent  data  inputs  and  interpolating  the  data  fcr  the  wave- 
length intervals  required  by  the  integration  subroutine.  IFLAG  can  take  on 
values  between  0 and  4 which  in  turn  keys  the  routine  to  read  the  plume 
intensity,  plume  transmission,  optical  transmission,  detectivity  and  solar 
irradiance  inputs,  respectively.  Once  this  data  has  been  read,  the  routine 
checks  to  see  whether  the  specified  wavelength  interval  for  the  passband  of 
interest  lies  outside  the  minimum  or  maximum  values  specified  by  the  inputs. 

If  the  wavelength  of  interest  lies  within  the  data  bounds  (as  it  should' 

YNEW  is  determined  using  a linear  interpolation  routine.  IFLAG  >5  allows 
for  interpolation  of  the  cloud/sky  radiance  tables.  Figure  27  shows  a 
logic  flow  for  this  subroutine. 

Subroutine  LOWEM.  This  subroutine  is  essentially  LOWTRAN  38  with 
path  emissions  added.  This  program  calculates  the  transmittance  and  ir- 
radiance of  the  atmosphere  from  350  inverse  centimeters  to  40,000  inverse 
centimeters  (0.25  to  28,5  microns)  at  30  inverse  centimeters  spectral 
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resolution  on  a linear  wavenumber  scale.  Refraction  and  ear;.!  curvature 
effects  are  included.  The  atmosphere  is  layered  in  1 kilometer  intervals 
between  0 and  25  kilometers,  5 kilometer  intervals  to  50  kilometers,  a 20 
kilometer  interval  to  70  kilometers  and  a 30  kilometer  interval  to  100 
kilometers.  There  are  six  model  atmospheres  available:  tropical,  mid-latitude 
summer,  mid-latitude  winter,  sub-arctic  summer,  sub-arctic  winter,  and  the 
1962  U.S.  Standard.  The  type  of  path  is  input  using  a code  from  1 to  6 to 
select  one  of  these  atmospheres.  It  is  also  possible  to  input  a new  atmos- 
phere such  as  radiosonde  data  in  which  case  a seventh  option  would  be  indicated. 
Aerosol  scattering  can  be  computed  and  there  is  a haze  model  available. 
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There  are  three  types  of  atmopsheric  paths  - vertical  or  slant 
paths  to  space,  vertical  or  slant  paths  between  two  altitudes,  and  a con- 
stant pressure  path.  The  output  of  this  subroutine  is  the  path  radiance  and 
atmospheric  transmission.  Subroutine  L0WEM  requires  Subroutine  POINT  which 
calculates  the  mean  refractive  index  above  and  below  a given  altitude  and 
interpolates  exponentially  to  determine  the  equivalent  absorber  amount  at 
that  altitude.  Also  required  for  L0WEM  is  Subroutine  ANGL  which  calculates 
the  initial  zenith  angle,  taking  into  account  refraction  effects. 

Subroutine  SCALR.  This  subroutine  adjusts  the  magnitude  of  the 
voltage  array  to  facilitate  the  evaluation  of  the  computer  printout.  The 
minimum  and  maximum  signal  levels  are  first  found.  If  the  data  is  complex, 
the  amplitude  is  used  for  finding  the  minimum  and  maximum.  The  minimum 
is  then  used  to  scale  the  data.  A multiple  of  10  is  calculated  such  that 
when  the  data  is  multiplied/divided  by  this  factor,  the  minimum  does  not 
exceed  100.  Using  the  factor  based  on  the  minimum,  the  voltages  are  dis- 
played. 


Since  the  minimum  is  used,  it  is  conceivable  that  the  maximum  value 
will  exceed  two  digits,  and  the  signal  level  will  be  displayed  with  an 
asterisk  (*).  If  this  is  the  case,  a new  scale  factor,  based  on  the  maximum, 
is  used  to  scale  the  voltage  array  again.  Glint  points  are  ignored  during 
these  scaling  calculations  and  are  always  displayed  with  asterisks. 


> 
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Subroutine  SKYCLD.  Subroutine  SKYCLD  is  one  of  the  more  complicated 
routines  within  IRSTEM.  It  is  here  that  the  signal  response  due  to  background 
sky  and  clouds  are  computed  and  substituted  into  the  array  (see  Fig.  28). 
There  are  six  types  of  sky  conditions  as  indicated  on  Fable  3.  IC1  indi- 
cates a clear  sky;  this  parameter  keys  the  program  to  compute  the  path  ra- 
diance only.  IC2  and  IC3  represent  hazy  conditions  with  23  and  5 km  visi- 
bility, respectively.  These  conditions  vary  from  the  clear  sky  computations 
in  that  solar  scattering  is  added  to  path  radiance.  The  last  three  types 
of  sky  conditions  are  clouds.  The  first  two  (100  m and  500  m clouds)  repre- 
sent the  case  when  the  observer  is  on  the  opposite  side  of  the  cloud  from 
the  sun  and  sees  partial  transmission  as  well  as  scattering  and  thermal 
radiation  from  the  cloud.  The  last  case  (10  km  cloud)  models  the  back 
reflection  from  an  optically  thick  cloud.  Table  3 indicates  the  dimensions 
of  the  solar  zenith  angle  (SZA)  solar  azimuth  angle  (AZA)  and  view  angle  (VA) 
which  must  be  known  to  properly  read  the  cloud  radiance  data  that  is  on  TAPE  8. 
The  specific  order  of  the  data  files  is  shown  in  the  right  hand  column  of 
Table  3. 

The  SKYCLD  computations  begin  by  calling  the  sun  location  subroutine 
to  compute  the  azimuth  and  elevation  angles  of  the  sun.  A series  of  anqles 
are  calculated  next  which  are  valid  for  all  of  the  sky  and  cloud  options  except 
10  km  (reflective)  cloud.  It  should  be  noted  that  the  cloud  normal  is  forced 
to  coincide  with  the  sun  angle  for  all  cases  except  the  10  km  cloud  to  simplify 
the  calculation  of  relative  angles  (see  Figure  7).  A series  of  flags  are  set 
before  preceding  into  the  signal  response  calculations  to  insure  the  proper 
sky  and  cloud  radiance  values  are  used.  A detail  flow  diagram  of  the  flag- 
setting algorithm  is  shown  in  Figure  29. 

The  initial  step  in  setting  the  flags  is  to  check  to  see  if  the 
clear  sky  conditions  prevail;  if  so,  the  next  few  steps  are  bypassed 
and  the  calculations  proceed  as  shown  in  Figure  28.  Next  a check  is  made  to 
see  if  the  10  km  thick  cloud  is  specified;  if  so,  the  number  of  unique  normals 
describing  the  cloud  are  read  followed  by  the  flags  and  data  limits  for  TAPE  8. 
If  neither  IC1  nor  IC6  are  applicable,  the  subroutine  DUMR  is  called  to  skip 
over  the  10  km  cloud  file,  and  a check  is  made  to  see  if  the  100  m cloud  is 
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now  under  consideration.  Again,  if  it  is,  the  flags  and  data  limits  are  set; 
otherwise,  the  dummy  read  is  called  again  to  skip  over  that  file  and  precede 
to  the  remaining  cloud  and  sky  types.  Control  is  transferred  from  this  part  of 
the  subroutine  after  all  6 types  of  sky  and  cloud  conditions  have  been  checked. 

Returning  to  Figure  28  we  see  that  the  major  DO  loop  for  the  number 
of  unique  normals  is  entered  following  the  flag  specifications.  The  number 
of  cells  of  the  background  array  containing  this  normal  are  read  and  a check 
is  made  to  see  if  the  10  km  cloud  is  being  analyzed  at  this  point.  If  it  is, 
the  elevation  and  azimuth  angles  defining  the  cloud  normal  are  read.  A check 
is  then  made  to  determine  if  the  clear  condition  is  being  analyzed;  if  so, 
the  next  several  steps  are  bypassed.  If  the  10  km  cloud  is  the  situation 
under  examination.  Subroutine  GETANG  is  called  to  compute  the  solar  angles 
in  terms  of  the  local  cloud  normal.  Then  after  some  more  preliminary  cal- 
culations, Subroutine  FIND  is  called  three  times  in  succession  to  identify 
the  radiance  data  files  which  most  closely  bound  the  specified  angles.  This 
is  necessary  to  insure  /hat  an  accurate  interpolation  of  the  radiance  data 
tape  is  made.  After  this  interpolation  is  completed,  the  next  major  DO  loop 
is  entered  to  compute  the  voltage  signals  arising  for  each  cloud  or  sky 
condition  within  the  background  scene.  LOWEM  is  called  and  the  voltage 
contribution  due  to  path  radiance  is  determined.  If  clear  sky  conditions 
are  being  analyzed,  these  are  the  only  contributions  to  the  background 
signal  level.  If  any  other  conditions  are  being  evaluated,  the  signal 
response  to  the  radiance  data  on  TAPE  8 is  computed  and  added  to  the  path 
radiance. 

The  only  remaining  SKYCLD  function  is  to  apply  the  computed  voltages 
to  the  appropriate  background  cells.  This  can  be  done  either  by  adding  the 
signals  to  an  area  specified  by  the  bounding  rows  and  columns  or  specifying 
the  location  of  each  sky  or  cloud  cell,  one  at  a time.  Control  is  then 
transferred  back  to  the  start  of  the  routine  to  process  the  next  sky  or 
cloud  type  until  all  6 types  have  been  processed. 
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The  sequence  of  READ  statements  within  SKYCLD  are  summarized  in  Table  4. 
The  first  time  through  the  subroutine,  checks  are  made  for  clear  sky.  If 
present,  the  number  of  cells  which  are  clear  and  the  minimum  and  maximum 
boundaries  or  the  individual  cell  locations  are  read.  On  the  next  pass  through 
the  subroutine,  a check  is  made  to  see  if  the  10  km  cloud  is  present.  If  it  is 
present,  the  number  of  calls  with  those  normals,  the  defining  angles  for  the 
cloud  normal,  and  the  cell  locations  are  read.  All  of  the  rest  of  haze  and 
thin  cloud  conditions  are  read  in  a manner  similar  to  clear  sky,  in  the 
sequence  shown  by  the  table. 


Subroutine  SUNLOC.  The  objective  of  SUN ICC  is  to  calculate  the  alti- 
tude (vertical  angle  measured  above  the  boric  true  bearing  (azimuth 

angle  measured  clockwise  from  true  north)  or  the  sun  as  a function  of  the 
geographical  position  of  the  observer  (latitude  and  longitude),  date  (month, 
day),  and  time  of  day  (Greenwich  Mean  Time).  To  develop  a solution,  one 
utilizes  the  celestial  triangle  as  reviewed  in  references  1 and  3.  The 
pertinent  equations  are  specified  in  Section  2.7. 

Subroutine  TGTSIG.  The  target  signar^.-q  -u:  • ie  computes  the  total 
integrated  intensity  for  each  passband  of  interest  and  sums  the  projected 
areas  of  the  overall  missile.  It  begins  by  computing  the  projected  areas  of 
the  tailpipe,  nose  and  missile  body.  Then  for  each  of  the  passbands, 
the  tailpipe,  missile  nose,  and  missile  body  signatures  are  computed  using 
the  blackbody.  The  plume  signature  then  is  computed  and  all  four  terms 
are  combined  to  give  the  overall  radiant  intensity.  The  equations  used 
in  this  routine  are  summarized  in  Section  2.7. 

Subroutine  TRAJR.  This  subroutine  locates  the  position  of  the 
target  on  the  background  with  respect  to  an  I,  J index  as  well  as  the  range 
from  the  sensor.  Currently,  this  trajectory  routine  assumes  that  the 
target  is  a point  source  flying  radially  inbound  with  an  aspect  angle  which 
does  not  change  with  time.  Since  these  computations  are  in  a separate  sub- 
routine, the  trajectory  information  can  be  fed  from  an  engagement  simulator 
if  desired. 


Table  4.  Read  Sequence  for  SKYCLD* 


NOCELS 

ICMIN/MAX,  JCMIN/MAX 
or 

NI(K) , NJ(K) 

NORM 

NOCELS 

THN , PHN 

ICMIN/MAX,  JCMIN/MAX 
or 

NI(K),  NJ(K) 

NOCELS 

ICMIN/MAX,  JCMIN/MAX 
or 

NI(K),  NJ(K) 


0.1  km  cloud,  0.5  km  cloud, 
23  km  haze  and  5 km  haze 
(in  turn) 


* Assuming  all  sky  and  cloud  types  are  present 
in  background  scene 
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3.3  Subroutine  Summary  Tables 

This  subsection  consists  of  tables  summarizing  the  inputs  and  outputs 
of  each  subroutine.  The  first  table  (Table  5)  is  the  cross-reference  map  of 
calls  to  each  subroutine.  Note  that  most  of  the  calls  come  from  either  the 
main  program,  IRSTEM,  or  the  SKYCLD  subroutine.  Tables  6-20  summarize  the 
key  parameters  within  each  subroutine.  Note  that  the  column  marked  "location" 
identifies  how  the  parameter  is  transferred  to  and  from  the  subroutine.  If  a 
common  block  is  used,  its  name  appears  in  capital  letters. 
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Table  5 . Subroutine  Call  Reference  Map 
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dy31NI 

Aagjna 


» 


Table  6 

Subroutine  FIND  Parameters 


Input  Parameter 

Location 

Description 

THO 

Call  Statement 

Desired  0O,  4>  or  0 on  successive 
calls  from  SKYCLD  (degrees). 

ICI 

Call  Statement 

Number  of  data  entries. 

SZA 

Call  Statement 

Value  of  0O.  $ or  9 for  which 
radiance  data  exists  (degrees). 

Output  Parameter 

Location 

Description 

J1 

Call  Statement 

Data  entry  index  for  SZA 
immediately  below  THO. 

J2 

Call  Statement 

Data  entry  index  for  SZA 
immediately  above  THO. 

D 

Call  Statement 

Relative  position  of  THO 

between  SZA  entries. 
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Table  7 

Subroutine  FIX  Parameters 


Input  Parameter 

Location 

Descripti on 

ISIGN 

Call  Statement 

ISIGN  equal s 1 for  an 
inverse  transform 

ISIGN  equals  -1  for  a 
ward  transform 

for- 

NUMBER 

Call  Statement 

NUMBER  equals  1 for  a 
plex  data;  0 for  real 

com- 

data 

U 

Call  Statement 

' 

Data  to  be  transformed 

Output.  Parameter 

Locati on 

Description 

U 

Call  Statement 

Transformed  data 
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Table  8 

Subroutine  FOL'RT  Parameters 


* 


1 

I 

I 


Input  Parameter 
DATA 

NN 

NDIM 

I SIGN 

I FORM 

WORK 

Output  Parameter 

DATA 


Location 
Call  Statement 

Call  Statement 

Call  Statemert 

Call  Statement 

Call  Statemert 

Call  Statement 

Location 


Descr ' t • on 

Data  to  be  transformed,  n-responds 
to  VARRY 

Dimension  array  describing  the  data 
dimensions,  i.e.,  NK'{1)  - 54,  NN(2)  ■ 
64,  sets  up  a two-dimensional  trans- 
form of  64  x 64  array 

The  dimensionality  of  the  transform, 
in  this  case,  2 

Integer  giving  direction  of  transform 
to  be  done  corresponds  to  ISIGN  of 
Subroutine  FIX 

Integer  parameter  describing  the  form 
of  the  data  which  corresponds  to  NUMBER 
of  Subroutine  FIX 

Dummy  variable  or  working  array  for  the 
fast  Fourier  transform 


Descript i on 


Formal  Parameter  Transformed  data  is  put  back  into  the 

input  data  array 


Table  9 

Subroutine  GETANG 


Input  Parameter 

Location 

Description 

ALT 

ANG 

Elevation  angle  of  sun  (rad). 

PHL 

GLT 

True  bearing  of  observer  as  seen 
from  the  cloud  (rad). 

PHN 

ANG 

True  bearing  of  cloud  normal  (rad). 

THL 

GLT 

Zenith  angle  of  observer  as  seen 
from  cloud  (rad) . 

THN 

ANG 

Zenith  angle  of  cloud  normal  (rad). 

TRUEB 

ANG 

True  bearing  of  sun  (rad). 

Output  Parameter 

Location 

Description 

TH 

ANG 

Angle  of  observer  from  cloud  normal 
(rad). 

THO 

ANG 

Angle  of  sun  from  cloud  normal  (rad 

PH 

ANG 

Angle  between  sun  and  observer  in 
plane  of  cloud  measured  from  reci- 
procal of  the  sun  direction  (rad). 

Subroutine  GLINT  Parameters 


Input  Parameter 

Location 

Description 

IN  IT 

GLT 

Flag  = 0 until  GLINT  is  called  then 
= 1. 

NOGLT 

Card  Read 

Number  of  GLINT  points. 

PSG 

Card  Read 

Probability  that  glint  does  not 
appear  in  a designated  cell  on  the 

scan. 

NG( 1,1) 

Card  Read 

Row  number  for  glint  I. 

NG( I ,2) 

Card  Read 

Column  number  for  glint  I. 

RGLINT ( I ) 

Card  Read 

Range  to  glint  I (m) . 

VEFF 

GLT 

Static  voltage  array  (v). 

Output  Parameter 

Location 

Description 

VEFF 

GLT 

Voltage  array  after  glint  is  super- 
imposed on  a specified  cell  (v). 

Table  11 

Subroutine  INTEG  Parameters 


Input  Parameter 
N 

Y(300) 

WAVEL  (100,3) 


Type/Location 


inscription 


INTP 


number  of  discrete  points  of  the 
function 


INTP 

CONST 


discrete  points  of  the  function 
to  be  integrated 

wavelength  integration  points 


Output  Parameter 

Y(l) 


Type/Location 

INTP 


Description  

result  of  the  integration. 
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Table  12 

Subroutine  INTER  Parameters 


» 


Input  Parameter 

Location 

Description 

I FLAG 

INTP 

flag  specifv’rig  whicn  spectral 
curve  is  being  read  and  interpolated. 

N 

INTP 

number  of  points  to  be  read  defining 
spectral  curve. 

NPW(3) 

INTG 

number  of  integration  points  in  each 
curve  for  each  given  waveband. 

WAVEL  (100,3) 

CONST 

wavelength  integration  points. 

XX  (200) 

Card  Read 

wavelengths  defining  the  spectral 
curve. 

YY  (200) 

Card  Read 

functional  values  at  given  wavelength 
values. 

Output  Parameter 

Location 

Description 

YNEW  (100,3) 

INTP 

interpolated  values  as  a function  of  _-vc 
inverse  centimeter  spacinq 

f 


12. 


Table  13 


Subroutine  LOWEM  Parameters 


Input  Parameter 

Type/Location 

Descri ption 

MODEL 

Call  Statement 

integer  to  select  the  model  atmosphere: 
model  1 - tropical 
model  2 - mid-latitude  summer 
model  3 - mid-latitude  winter 
model  4 - sub-arctic  summer 
model  5 - sub-arctic  winter 
model  6 - 1962  U.S.  Standard 
model  0 - horizontal  path  when  meteoro- 
logical data  are  used 
model  7 - new  model  atmosphere 

I HAZE 

Call  Statement 

integer  to  indicate  aerosol  attenuation 
IHAZE  0 - no  aerosol  attenuation 

IHAZE  1 - aerosol  attenuation  calculated. 

ITYPE 

Call  Statement 

integer  to  indicate  the  type  of  atmos- 
pheric path: 

ITYPE  3 - vertical  or  slant  path  to  space 
ITYPE  2 - vertical  or  slant  path  between 
two  altitudes 

ITYPE  1 - horizontal  (constant  pressure) 
path 

LEN 

Call  Statement 

LEN  1 - for  a downward  looking  when  two 
paths  are  possible  (normally 
set  to  0) 

JP 

Call  Statement 

print  option 

IM 

Call  Statement 

IM  1 - for  radiosonde  data  (normally  set 
to  0) 

Ml,  M2,  M3 

Call  Statement 

these  integer  values  are  used  to  change 
the  temperature,  the  water  and  the  ozone 
altitude  profiles  (normally  set  to  0) 

ML 

Call  Statement 

ML  1 - for  radiosonde  (normally  set 
to  0) 

IEMIS 

Call  Statement 

integer  to  indicate  mode  of  calculations: 
IEMIS  0 - transmission  mode 

IEMIS  1 - radiance  mode.  (Always  in  IRSTEM.) 
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Subroutine  LOWEM  Parameters 


Input  Parameter 

Type/Location 

Description 

» 

RO 

Call  Statement 

range  (km) 

HI 

Call  Statement 

height  of  observer  in  kilometers 

H2 

Call  Statement 

source  altitude  in  kilometers 

i 

ANGLE 

Call  Statement 

zenith  angle  at  HI  in  degrees 

RANGE 

Call  Statement 

pathlength  in  kilometers 

BETA 

Call  Statement 

earth  center  angle 

» 

VIS 

Call  Statement 

visual  range  at  sea-level  in  kilometers 

WAV  El 

Call  Statement 

lower  bound  of  calculations  in  microns 

1 

WAVE2 

Call  Statement 

upDer  bound  of  calculations  in  microns 

DELAM 

Call  Statement 

increment  of  calculations  in  microns 

IXY 

Call  Statement 

ISP 

Call  Statement 

integer  to  indicate  previous  calculations 

1 * 

have  bee:’  June: 

ISP  - need  first  call  to  Subroutine 

l:x 

ISP  1 - lc  c 're  than  one  call  to 

LOWEM 

N 

I WRITE 

Call  Statement 

print  FLAG 

UK 

Call  Statement 

not  used 

I STATE 

Call  Statement 

not  used 

<11 

IAFLG 

Call  Statement 

not  used 

TOTRAN 

Call  Statement 

array  names  for  transmission  values 

LW 

Call  Statement 

dimension  counter 

1 

NSAI 

Call  Statement 

0 - Burch  H:0  model,  1 - White  model 

Output  Parameter 

Type/Location 

Description 

TAUA 

LOW 

transmissivity 

• 

APRA 

LOW 

path  radiance 
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Table  14 

Subroutine  SCALR 


Input  Parameter 

Location 

Description 

EFF  (64,  64,  1) 

Call  Statement 

Signal  response  array  (v) 

I MAX 

PRT 

Number  of  rows  in  scene. 

JMAX 

PRT 

Number  of  columns  in  scene 

L 

Call  Statement 

Waveband 

LAVEL 

Call  Statement 

Label  for  array  printout 

NG  (40,  2) 

PRT! 

Location  of  glint  points 

NOGLT 

PRT1 

Number  of  glints 

Output  Parameter 
IBK6  (64,  64) 


Location  Description 

PRT  Scaled  background  array. 


Table 

15 

Subroutine 

SKYCLD 

Input  Parameter 

Location 

Description 

ALT 

ANG 

Elevation  angle  of  sun  (rad). 

APR 

LOW 

Path  radiance  (mw/sr/cnr/um) . 

AZA 

Data  Statement 

Azimuth  angle  between  sun  and  obser- 
ver in  cloud  plane  for  which  radiance 
values  exist  (rad) . 

D 

MISC 

Optics  diameter  (m) 

DI,  DJ,  DK 

Subroutine  FIND 

Relative  position  of  THO,  PH  and 

TH  between  the  bounding  data  entries. 

DU.M  (300)  or 

YNEW  (300 

IilTP 

Interpolated  integrands  for  signal 
response  calculations. 

0V(  3) 

INTG 

Wavelength  interval  (um) 

ELVS 

ITRAJ 

Zenith  angle  of  sensor  (rad). 

11,  12 

Su  ' 

Data  entry  indices  for  SZA  immediately 
below  and  above  THO. 

IC1-IC6 

GLT 

Flags  indicating  the  sky  or  cloud 
condi tion. 

ICMIN,  ICMAX 

Card  read 

Rows  bounding  a given  sky  or  cloud 
condition. 

Jl,  J2 

Subroutine  FIND 

Data  entry  indices  for  AZA  immediately 
below  and  above  PH. 

JCMIN,  JCMAX 

Card  read 

Columns  bounding  a given  sky  or 
cloud  condition. 

Kl,  K2 

Subroutine  FIND 

Data  entry  indices  for  VA  immediately 
below  and  above  TH. 

NBANDS 

INTG 

Number  of  passbands. 

NI(K),  NJ(K) 

Card  Read 

t h 

Row  and  column  location  of  the  Ku 
cell  of  a given  sky  or  cloud  condition 

NOCELS 

Card  read 

Number  of  background  cells  with  a 
given  normal. 

±J 


Table  15  (Continued) 


Input  Parameter 

Location 

Description 

NORM 

Card  Read 

Number  of  unique  normals  for  a qiven 
sky  condition  (set  to  1 internally 
except  for  10  km  clouds). 

PH 

Internal  calculation 

Azimuth  angle  between  sun  and  observer 
in  cloud  plane  (rad)  - from  ANG  only 
for  10  km  cloud. 

PHL 

GLT 

True  bearing  of  1 ine-of-sight  (rad). 

PHN 

Card  read 

Azimuth  angle  of  cloud  normal  (rad). 

RB(40) 

IMPT 

Range  to  cloud  (m) 

RD 

Responsivity  (v/w) 

SZA 

Data  statement 

Solar  zenith  angle  wrt  cloud  normal 
for  which  radiance  data  exists  (rad). 

TAUA 

LOW 

Atmospheric  transmissivity. 

TAUO 

GLT 

Optical  Train  Transmissivity. 

TH 

Internal  calculation 
or  ANG 

View  angle  of  observer  wrt  cloud 
normal  (rad)  from  ANG  only  for  10  km 
cloud. 

THN 

Card  read 

Zenith  angle  of  cloud  normal  (rad). 

THO 

Internal  calculation 
or  ANG 

Solar  zenith  angle  wrt  cloud  normal 
(rad)  from  ANG  only  for  10  km  cloud. 

TRUEB 

ANG 

True  bearing  of  sun  (rad). 

VA 

Data  statement 

Azimuth  angle  between  sun  and  observer 
in  cloud  plane  for  which  radiance  values 
exist  (rad). 

XLAML(L) ,XLAM(U) 

GLT 

Upper  and  lower  wavelengths  for  each 
waveband  (urn). 

YR 

Data  statement  or 

TAPE  8 

Wavelengths  for  which  radiance  data 
exists  on  TAPE  8 (pm)  or  the  radiance 
values  from  the  tape  (mw/sr/enr/ym) 

73 


Subroutine  SUNIOC 


Input  Parameter 

Location 

Descri ption 

MON 

GLT 

Month  oT  year  (1  - 12) 

DAY 

GLT 

Day  of  month  (1  - 31 ) 

HR 

GLT 

Time  of  day  (GMT) . 

RLAT 

GLT 

Latitude  of  interest  (degrees)  - 
North  is  +,  South  is  - . 

RLONG 

GLT 

Longitude  of  interest  (degrees) 

Output  Parameter 

Location 

Description 

ALT 

ANG 

Altitude  of  sun  from  horizon  (rad) 

TRUEB 

ANG 

True  bearing  of  sun  clockwise  from 
True  North  (rad). 

I 


I 


Table  17 

Subroutine  TGTSIG 


Input  Parameter 

Location 

Description 

ASPECT 

ITRAJ 

View  angle  of  target  wrt  its  center 
line  (degrees) 

DP 

MISS 

Tailpipe  diameter' (m) . 

DT 

ITRAJ 

Target  diameter  (m). 

ETE 

MISS 

Tailpipe  emissivity. 

ETN 

MISS 

Nose  emissivity. 

ETS 

MISS 

Skin 'emissivity. 

TE 

MISS 

Tailpipe  temperature  (°K) 

TL 

ITRAJ 

Target  length  (m). 

TN 

MISS 

Nose  temperature  (°K). 

TS 

MISS 

Skin  temperature. 

TOP  (100,  3) 

IMPT 

Plume  spectral  intensity  (w/m/pm) 

TNL 

MISS 

Length  of  missile  nose  (m) 

TGTAU  (100,  3) 

IMPT 

Plume  spectral  transmissivity. 

WAVEL  (100,  3) 

CONST 

Wavelength  of  interest  (urn) 

Ouptut  Parameter 

Location 

Description 

TJ 

ITRAJ 

Total  target  intensity  (w/sr^um) 

AT 

ITRAJ 

Total  target  cross-section  (m2) 

oO 


Table  18 


Subroutine  TRAJR  Parameters 


Input  Parameter 

ELVS 

Type/ Location 

ITRAJ 

IMAX 

PRT 

JMAX 

PRT 

PHITS 

ITRAJ 

RS 

ITRAJ 

T 

ITRAJ 

TF 

ITRAJ 

THTATS 

ITRAJ 

VT 

ITRAJ 

Outout  Parameter 

Type/Location 

II 

ITRAJ 

JJ 

ITRAJ 

R 

ITRAJ 

Description 

theta  position  of  the  sensor  (rad) 

number  of  rows. 

number  of  columns. 

azimuth  position  of  the  target  (rad) 
initial  range  of  the  target  (m) 
engagement  clock  time  (sec) 
frane  time  (sec) 

theta  position  of  the  target  (rad), 
velocity  of  the  target  (m/sec). 

Description  _ _ 

I index  position  of  Ui scl. 

J index  position  of  the  target 
range  of  the  target. 


▲ 


Function  BLKBDY 


Input  Parameter 

Location 

Description 

Cl 

Call  Statement 

Constant  * 3.713  x 104  (w  u4cm’2 

C2 

Call  Statement 

Constant  = 1 .4388  x 104(u°K) 

TT 

Call  Statement 

Temperature  (°K) 

XL 

Call  Statement 

Wavelength  (y) 

Output  Parameter 

Location 

Description 

BLKBDY 

Call  Statement 

Blackbody  emittance  (w  cm"- 

Table  20 
Function  INTERP 


Input  Parameter 
ARG 

N 


Location 
Call  Statement 


Description 

Aspect  angle  (degrees). 


Call  Statement  Number  of  points  defining  plume 

intensity  curve  (nominally  - 12) 


X 


Call  Statement  Function  describing  the  curve. 


Output  Parameter 

Location 

Description 

Y 

Call  Statement 

Function ' 
angle  of 

' s value  at  the  aspect 
interest. 

Section  4 

USER'S  INSTRUCTIONS 

In  this  section  code  design  information  is  provided  to  assist  the 
IRSTEM  user.  The  first  part  of  this  section  reviews  the  input  requirements 
and  the  output  options.  The  key  parameters  used  in  the  code  are  summarized 
on  a table  which  includes  the  parameter,  its  symbol  in  the  text  and  where 
it  is  defined,  stored  and  used.  This  table  combined  with  the  List  of 
Symbols  and  Subroutine  Tables  in  Section  3 give  a complete  cross-reference 
for  each  of  the  major  parameters.  Finally,  the  section  is  completed  with 
a sample  program  and  test  case. 

4. 1 Inputs/Outputs 

There  are  36  sets  of  inputs  required  to  operate  IRSTEM,  33  formal 
inputs  and  three  namelists.  These  inputs  include  a combination  of  tapes, 
formal  input  cards  and  name lists  structured  to  provide  maximum  flexibility 
for  the  user.  Some  of  the  inputs  are  simply  a single  parameter  while  others 
are  full  arrays  of  many  thousands  of  elements.  There  are  no  default  options 
in  IRSTEM;  so  the  user  must  provide  data  for  each  of  these  inputs  whenever 
the  program  is  exercised. 

Table  21  summarizes  the  input  stream.  All  of  the  inputs  are  assumed 
to  be  cards  or  files  of  data  in  card  image  format  unless  otherwise  noted  as 
tapes  or  namelists.  Each  entry  shows  the  requisite  format,  the  calling 
routine  (if  other  than  the  main  program),  the  parameter's  units  and  a des- 
cription of  each  parameter  including  the  flag  options.  The  paragraphs  below 
discuss  these  inputs,  in  turn,  and  provide  additional  information  for  the  user 
beyond  what  is  in  the  table. 

Program  operation  begins  by  specifying  the  background  array  size  and 
whether  the  debug  option  for  additional  printout  is  desired.  It  should  be 
noted  that  in  addition  to  specifying  IMAX  and  JMAX,  arrays  other  than  64  by 
64  require  the  user  to  alter  several  dimension  statements  and  common  blocks. 
The  required  modifications  are  shown  in  Table  22.  Most  of  the  changes  occur 
within  the  main  program;  however,  five  of  the  subroutines  do  require  array 
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Table  21 
IRSTEM  Inputs 


Input  1 - 315 


I MAX 

- 

Number  of  rows  in  background  array. 

JMAX 

- 

Number  of  columns  in  background  array. 

NBUG 

- 

0 -*•  standard  output 

Input  2 - 5A6.I5 

(Tape  1) 

1 -*■  expanded  output 

HOL 

- 

Title 

MXCODE 

- 

Maximum  number  of  different  codes  on  the  back- 

Input 3 - 1515  (Tape  1) 

ground  tape. 

IBKG(I.J) 

- 

Array  of  background  codes  from  1 to  IMAX 

Input  4 - I5.4F10 

.3 

and  1 to  JMAX. 

KEY(I) 

- 

Integer  value  corresponding  to  background  code. 

T8(I) 

°K 

Temperature  for  Ith  background  code. 

EB(I) 

- 

Emissivity  for  Ith  background  code. 

RHO(I) 

- 

Diffuse  reflectivity  for  Ith  background  code. 

RB(I) 

m 

Range  to  Ith  background  code. 

Namelist  SENSOR 

XLAML(I) 

pm 

Lower  bound  of  the  wavebands  of  interest  [ 

Upper  bound  of  the  wavebands  of  interest  ) 

XLAMU(I) 

pm 

TF 

sec 

Scan  time 

D 

m 

Diameter  of  the  optics 

DX 

m 

Width  of  the  detector 

DY 

m 

Length  of  the  detector 

OVL 

m 

Amount  of  overlap  between  detectors. 

FL 

m 

Focal  length 

FO 

hz 

Cutoff  frequency 

ELVS 

degrees 

Elevation  angle  of  the  background 

VNOISE 

V 

RMS  noise  level 

RDPK 

v/w 

Peak  responsivity 

DSTRPK  cm 

hzji  w"1 

D*  or  detectivity 

2,3 
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Table  21  (Continued) 
IRSTEM  Inputs 


Input  5 - 215,  2E10.3,  A8 


MODET 

NBANOS 

XLAML(l) 

XLAMU(l) 

NAMET 


un 


Type  of  trajectory,  trajectory  inode  (not  used) 
Number  of  wavebands  to  be  considered  (1,  2 or  3) 
Lower  bound  of  first  waveband 
Upper  bound  of  first  waveband 
Name  of  the  target. 


Input  6 - 415,  4F10.3,  15 
MODEL 
I HAZE 

ITYPE 
I EMI  S 

HS  in 


Integer  code  for  type  of  atmosphere  (see  Table  13) 

Integer  to  indicate  if  aerosol  scattering 

is  to  be  used  (0-no  aerosol,  l-^aerosols  included) 

Type  of  atmospheric  path  (see  Table  13) 

Mode  code  to  indicate  transmittance  or 

radiance  (must  be  1 if  path  radiance  is  to  be  included) 

Observer  altitude 


DV(  I ) 
NSAI 


■m  Integration  intervals  (I=1,NBANDS) 

NSA 1=1  indicates  that  White  water  vapor 
continuum  is  to  be  used  in  place  of  the 
Burch  model. 


Input  7 - 15  from  INTER 

N - Number  of  (A,  Jp)  pairs 

Input  8 - 4(F12.3,L8.2)  from  INTER 

A(I)  urn  Wavelength 

TJP(I)  w/sr/-,.m  Target  plume  intensity 


Input  9 - 15  from  INTER 
N 


Number  of  (A, 


r ) Pairs 


Input  10  - 9XF7.4,  52XF6.4 
A(I)  um 

TGTAU(I) 


from  INTER 
Wavelength 
Plume  transmission 


1 = 1 ,N 


Input  11  - 15  from  INTER 
N 


Number  of  (A,  r^)  pairs) 


I 


Table  21  (CoMf  tnued) 
IRSTEM  Inputs 


Input  12  - 3F10.J 

from  INTER 

MU 

m 

Wavelength  | , . 

} i \ » ii 

Optical  Transmission  > 

TAU0( I ) 

- 

Input  13  - 15  from 

INTER 

N 

- 

Number  of  ( ' , D* ) pa  i rs 

Input  14  - 3F10.3 

from  INTER 

Ml) 

1/2  -1 
z W 

Wavelenoth  | , , ,, 

> i • i , r< 

D*  ( I ) cm  h 

Detectivity  » 

.nput  15  - 15  from 

INTER 

N 

- 

Number  of  ( \,  M ) pairs 

Input  16  - 3F10.3 

from  INTER 

Ml) 

Wavelength  | 

sr(h 

w / car 

Solar  irradianco  it  oa  to’s  surface  1 

Namelist  TRAD 

RS 

m 

Initial  ramie  ot  target 

TV 

m/sec 

Velocity  ot  target 

DT 

m 

Diameter  ot  target 

TL 

m 

Length  of  target  bodv 

THTAIS 

degrees 

Elevation  of  tamet 

PH  ITS 

degrees 

Azimuth  of  target 

TNL 

m 

Length  of  Nose 

TPL 

m 

Plume  Length  (not  used  at  this  time) 

OP 

m 

Diameter  of  the  exhaust  pipe 

TE 

°K 

Temperature  of  exhaust 

TN 

°K 

Temperature  of  the  nose 

TS 

°K 

Temperature  of  the  bodv  of  the  missile 

Table  21  (Continued) 
IRSTEM  Inputs 


Namel i st 

TRAJ 

(Continued) 

ETN 

- 

Emissivity  of  the  nose 

ETS 

- 

Emissivity  of  the  body  of  the  target 

ETE 

- 

Emissivity  of  the  exhaust 

ASPECT 

degrees 

Aspect  angle  of  the  target  relative  to  the 
line  of  sight 

TMAX 

sec 

Maximum  time  of  the  encounter 

Namel ist 

SUN 

RLAT 

degrees 

Latitude  of  the  observer 

RLONG 

degrees 

Longitude  of  the  observer 

MON 

- 

Month  (1-12) 

OAY 

- 

Day  (1-31) 

HR 

- 

Hour  of  the  day  (GMT) 

PHL 

degrees 

Locates  the  observer  with  respect  to  true 

North  in  a counterclockwise  direction  as 
seen  from  the  target 

Input  17 

- 15 

from  SKYCLD 

NOCELS 

- 

Number  of  cells  with  clear  sky  (Flag  for  Input 
If  N0CELS=0,  block  data  is  read) 

Input  18 

-1615  from  SKYCLD 

ICMTN 

j 

1 

| Row  limits  for  clear  sn>- 

ICMAX 

] 

JCMIN 

( 

f Column  limits  for  clear  skv 

JCMAX 

1 

) 

or 

Input  13 

- 1615  frcm  SKYCLD 

NI(K) 

K = 

1.  NOCELS 

Cell  by  cell  specification  of  clear  skv 

NJ(K) 

) 

positions 

Input  19 

- 15 

from  SKYCLD 

NORM 

- 

Number  of  unique  normals  for  10  km  thick  cloud 

as 


Table  21  (Continued) 
IRSTEM  Inputs 


Input  20  - 15  froiit  SKYCIO 
NOCELS 


Number  of  cells  with  10  km  thick  cloud. 


Input  21  - 2E10.3  from  SKYCLD 

THN  degrees  Zenith  angle  of  cloud  normal 

PHN  degrees  Azimuth  angle  of  cloud  normal 

Input  22  - 1615  from  SKYCLD 

Same  as  Input  13  for  10  km  thick  cloud 


Input  23 

Input  24 

Input  25 
Input  26 

Input  27 
Input  23 


Same  as  Inputs  17  and  18  for  0.1  km  thick  cloud 


| Same  as  Inputs  17  and  18  for  0.5  km  thick  cloud 

( Same  as  Inp.its  17  and  13  for  '3  km  haze 


Input  29 
Input  30 


Same  as  Inputs  17  and  18  for  5 km  haze 


N0GLT 

PSG 


15,  El 0.3  from  GLINT 

- Number  of  cells  with  potential  glints 

Probability  that  a glint  will  not  occur  on 
any  given  scan. 


Input  32  - 1615  from  GLINT 


NG(IJ  1 = 1,  N0GLT 
NG( I ,2) 


Row  position  of  Ith  glint 

k L. 

Column  position  of  Pn  glint 


Input  33  - 3E10.3  from  GLINT 

RGLNT(I)  m Range  to  each  glint  point 

1=1,  NOGLT 


Table  22 

Changes  Required  if  IMAX  or  JMAX  Do  Not  Equal  64 


modi fications  before  the  program  will  function  properly.  To  make  the  modi- 
fications, the  user  changes  the  array  dimension  from  64  to  tne  appropriate 
number.  It  is  not  required  that  this  number  be  a multiple  of  2,  but  the  Fast 
Fourier  Transform  will  ODerate  more  efficient! v if  it  is  a multiple  of  2. 

Inputs  2,  3 and  4 define  the  background.  After  the  header  information 
is  provided,  the  background  codes  (one  for  each  cell'  are  read  as  Input  3. 

Then  key  cards  are  read  (one  for  each  background  type'  to  associate  each  code 
with  the  desired  temperature,  emissivity,  reflectivity  and  range. 

The  next  3 sets  of  inputs  provide  information  about  the  sensor  and 
the  atmospheric  conditions  desired  for  the  engagement.  These  inputs  are 
fixed  parameters  that  do  not  vary  with  wavelength.  Namelist  Sensor  provides 
information  on  the  hardware  design  and  Input  6 is  a shortened  10WTRAN  input 
list  required  by  LOWEM.  The  rest  of  the  LOWEM  parameters  are  shown  in  the 
subroutine  table  (Table  13).  These  parameters  are  set  to  default  values 
internally  in  the  program. 

Inputs  7 through  16  are  the  spectrally  dependent  parameters  which 
must  be  read  in  the  order  shown.  Two  sets  of  inputs  are  needed  for  each 
parameter.  The  first  card  specifies  the  number  of  data  points  and  the 
re-ai.'.'-g  cards  give  the  wavelength  and  cor  ■ >:  ending  value  of  the  functional 
for  eav.;i  data  point.  These  inputs  are  called  from  Subroutine  INTER  which 
then  interpolates  the  function  as  required  for  the  integration  subroutine. 

Two  additional  namelists  follow  Input  16.  The  first  is  used  to  describe 
the  engagement  dynamics  and  input:  ''squired  for  ‘he  target  signature  cal- 
culation (Namelist  TRAJ).  The  second  specifies  the  location  of  the  engage- 
ment with  parameters  necessary  to  compute  the  sun’s  location  (Namelist  SUN). 

Inputs  17  through  30  are  associated  with  the  cloud  radiance  model. 

The  cloud  inputs  must,  of  course,  be  correlated  with  the  background  data 
read  at  the  start  of  the  program;  however,  these  cards  give  the  user  the 
flexibility  of  creating  any  type  of  cloud  or  clear/haze  air  conditions  he 
wishes.  The  required  read  sequence  is  clear  sky,  10  km  thick  clouds,  0.1  km 
thick  clouds,  0.5  km  thick  clouds,  23  km  hate  and  finally  5 km  haze  if 
present  in  the  background  description.  In  general,  there  are  two  inputs 
required  for  each  sky/cloud  condition.  The  first  is  the  number  of  cells 
whose  location  will  be  individually  specified  as  containing  this  sky  or 


cloud  condition.  The  next  input  then  specifies  the  location  of  these  cells. 

If  the  sky/cloud  condition  can  be  conveniently  described  as  a rectangular 
section,  the  user  has  the  option  of  specifying  the  bounds  of  the  rectangle 
instead  of  specifying  the  location  of  individual  cells.  However,  the  user 
cannot  specify  a combination  of  area  and  individual  cell  locations.  The 
above  read  sequences  are  valid  for  all  the  sky  and  cloud  conditions  with  the 
exception  of  the  10  km,  reflective  cloud.  Here,  the  first  input  is  the 
number  of  unique  cloud  normals  used  to  describe  the  cloud.  Then  the  number 
of  cells  with  those  normals  are  given  and  defined  by  their  zenith  and  azi- 
muth angles.  Finally,  the  location  of  these  cells  are  read  in  the  same 
manner  as  before  (i.e.,  by  block  areas  or  individual  cell  locations).  In- 
puts 20  through  22  are  repeated  then  for  each  of  the  normals  specified  in 
Input  19. 

The  final  three  inputs  are  associated  with  the  sun  glint  computations. 
Input  31  identifies  the  number  of  cells  containing  potential  sun  glints  and 
the  probability  with  which  the  sun  glints  will  occur.  Note  that  this  prob- 
ability is  applied  equally  to  all  glints.  Next,  each  glint  location  must  be 
specified  on  a cell -by-cell  basis  (there  is  no  option  for  specifying  sun 
glint  areas  at  this  time).  Finally  the  range  associated  with  each  of  the 
sun  glints  is  specified.  Normally  this  range  will  correspond  to  the  back- 
ground range  given  in  Input  4. 

The  IRSTEM  outputs  are  summarized  in  Table  23  which  shows  the  para- 
meter in  the  order  of  its  appearance.  Parameters  which  are  preceded  with  an 
asterisk  indicate  debug  options  which  are  printed  only  if  NBUG  > 0. 

The  output  sequence  generally  follows  the  input  sequence  with  the 
background  conditions  specified  first  followed  by  descriptions  of  the  sensor, 
target  and  environmental  conditions.  Note  that  when  Namelist  SENSOR  is 
printed,  the  wavelength  values  for  the  primary  band  are  not  repeated  since 
they  are  shown  at  the  beginning  of  the  output  listing.  The  outputs  which  are 
printed  on  each  iteration  of  the  time  loop  are  (1)  the  total  projected  area 
at  that  time,  (2)  the  target's  range  and  (3)  the  response  voltage  of  the 
target  cell  for  each  waveband.  The  signal  response  from  each  detector, 

Vn ( t ) , is  recorded  on  Tape  12  for  later  signal  processing  analysis  but  is 
not  printed.  To  more  fully  understand  the  output  sequence  and  format,  the 
reader  should  examine  the  sample  output  given  at  the  end  of  this  section. 
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Parameter 
I MAX,  JMAX 

HOL 

KEY(I) 

TB(I) 

EB(I) 

RHO(I) 

RB(I) 

*!8KG( I ,J) 


Uni  ts 


Input  5 . 

XLAML(J) 

f J=2 , uin 

XLAMU(J) 

( riBANDS 

Input  6 

0 

1 /2 

cm  hz ' ~ 

V') 

v/w 

H ( '• ) 

l:  Cl 

a 

Name  1 j s t 

i 

Name! i st 

SPIN 

-1 


*IBKG( I , J) 
■*VE-C'  I ,J  ,L 

^Tp(  ) 

JM  (\) 

Js  (\) 

Jp  (X) 

JT  (\) 


Table  23 
IRSTEM  Outputs 

Description 

Number  of  rows  and  columns  in  background  array. 
Title 


Descriptors  for  each  background  code. 


Input  oackground  array 

Tarnet  information  and  passbands 

LOWEM  inputs 


Spectrally  dependent  inputs. 


Background  array  with  sequenced  numbers. 

Jlgnal  response  to  static  :c..  .. .glints, 


w/sr,  . in  Target  signature  including  component  parts 


Table  2 3 (Continued) 
IRSTEM  Outputs 


? a rameters 
AN 

AS 

A_ 

» 

L 

A 

I 


rVEFF( I ,J,L) 


7VEFF 


VEFF  l 


•EFF 


Units 


sec 

_2 


m 

v 


Description 


Projected  target  areas. 


Scaled  scene  response. 
Scaled  transformed  array. 


Target  status  presented  for  each  iteration  of 
the  time  loop. 


Scaled  scene  array  after  the  sensor  distortions 
are  added. 

Final  signal  response  for  each  of  the  detectors 
is  written  on  TAPE  12. 
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4 . 2 Parameter  List 

A List  of  Symbols  is  given  in  the  front  of  this  report  identifying 
* the  symbols  used  in  the  theory  section.  In  Section  3.2  a series  of  tables 

(Tables  6 through  20)  are  provided  summarizing  the  inputs  and  outputs  associ- 
ated with  each  of  the  subroutines  and  in  Section  4.1  above,  the  inputs  and 
outputs  required  to  operate  the  code  are  identified.  In  this  subsection 
9 the  cross-references  for  the  key  parameters  is  completed  by  identifying  each 

of  the  important  code  symbols  as  they  appear  in  IRSTEM  in  alphabetical 
order. 

9 Each  of  these  key  parameters  are  listed  on  Table  24  which  includes 

the  corresponding  report  symbol,  its  units,  where  the  symbol  is  defined,  its 
common  block,  and  which  are  the  primary  routines  in  IRSTEM  in  which  it  is 
used.  Under  the  "Defined"  column  there  are  four  options.  Words  which  appear 
9 in  all  caps  (such  as  SUNLOC)  indicate  the  parameter  is  defined  in  a subroutine. 

If  the  word  begins  with  a capital  letter  (such  as  Sensor),  this  indicates 
a namelist  read.  Formal  read  inputs  are  indicated  as  Input  x and  are  based 
on  the  read  sequence  shown  in  Table  21.  Finally,  some  of  the  parameters 
9 are  defined  in  data  blocks.  The  dimensions  of  the  variables  are  qiven  in 

the  "Description"  column  where  apolicable,  and  if  the  array  is  a complex  one, 
this  also  is  noted. 


r 


35 


Table  24  (Continued) 


TRA,1R  ITPAJ  j IRSTEM  Row  in  which  tarqet  is  located. 


IRSTEM  INTP  INTER  Number  of  points  to  be  inter- 

polated. 


TAPE8  - - SKYCLD  Cloud  radiance  data. 
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4 . 3 Sample  Frog ram 

In  this  section,  a sample  program  is  presented  including  *he 
computer-generated  1’st  of  the  input  and  output  stream.  The  us 

test  case  is  designed  to  illustrate  the  use  of  the  parameter  •••rained 
in  the  input  output  tables  (Tables  21  and  23)  and  is  not  designed  to 
test  all  of  the  code's  features . 

Table  25  shows  the  input  stream  associated  with  this  sample 
program.  The  target  is  flying  radially  inbound  toward  a < • no ’o-col or 
sensor.  The  background  is  characterised  by  three  unique  elements;  a 
thick  cloud,  a 295°k  graybody  at  50  km,  and  a 300*" K graybody  at  20  km. 
It  should  be  noted  that  trie  allocation  of  these  background  elements  to 
the  4096  element  array  is  not  shown  on  this  input  list  since  if  is 
read  from  a separate  source  (Tape  1'.  It  also  should  he  noted  that 
the  standard  plume  transmission  tabulated  input  (see  Figure  l1'' 
of  172  points  is  not  shown  here.  Instead,  a few  of  the  beginning  and 
end  points  are  shown  along  with  the  ^ values  which  bound  the  3.S  - 
4.2..m  passband  selected  for  this  test  case.  Similar!;,  . only  toe  first 
few  lines  of  the  solar  irradiance  spectrum  (Figure  • ar->  snown  in  he 
interest  of  brevity.  Also  not  shown  are  a large  number  of  inputs  re- 
quired by  the  propagation  model.  These  inputs  are  the  standard  LOW’vAN 
inputs  covered  in  reference  4.  The  last  few  lines  of  the  input  stream 
define  the  cloud  angles,  position,  and  the  two  glint  points. 

In  Table  26  the  corresponding  outputs  for  ttr  . test  case 
are  shown  as  they  appear  for  the  NBUG-0  with  two  except'  s.  ■ ■'•>:. 
the  time  loop  is  truncated  after  three  iterations  (t-0,  1 and  2 
seconds)  even  though  TMAX*10  seconds.  Second,  several  arrav  print- 
outs are  included  which  normally  are  printed  only  if  NBUG-0.  These 
additional  arrav  plots  are  (1)  the  background  inputs  after  index 
normal isation,  (2  and  3)  the  signal  response  after  the  target 
has  been  added,  for  the  minimum  and  maximum  array  normalisation 
but  before  the  glint  and  system  noise  are  included  and  (4'  one  of 
the  transformed  arrays.  Many  more  arrays  are  printed  if  NBUG'O  out 

sfiow  i v".  '•  urse.  these  arrays  contain  no  re  inr  crma- 
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C i on  than  is  required  since  only  those  rows  corresponding  to 
detector  centers  are  stored  on  Tape  12  and  subsequently  used  in 

accessing  ana1 /set.  With  respec*  to  the  "'.tindard 
e - . so  c sr’ntouts,  it  arcald  -c  noted  tnat  tne  target's 
urojected  trea  and  its  soectral  intensity  are  repeated  even  thougn 
tney  do  no;  change  for  tins  radial ly-inbound  trajectory.  These  print- 
outs nave  been  included  in  anticipation  of  time-varying  signatures 
and  trajectories  to  be  incoroorated  in  later  versions  of  IRSTEM. 


Table  25 

Sample  Input  Stream 
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Tar-le  25  (continued) 
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Table  25  (continued) 
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Taole  26  (continued) 


Cl  RANGEIfll 
HAItSFCRnE}  VAll'E 
1 : 3 4 5 4 

1 5 5 J 3 5 5 

2 5 3 5 5 5 5 

3 5 5 3 3 5 5 

4 5 5 5 5 3 3 


11  3 3 3 

12  5 5 5 

13  5 5 3 

14  5 5 3 

13  5 3 5 

14  5 5 3 

17  5 5 5 

18  5 3 5 

H 5 3 3 

20  5 3 5 

21  3 5 5 

22  5 5 3 

23  5 5 3 

24  3 5 3 

23  5 3 3 

24  5 3 3 

27  3 5 5 

28  3 3 5 

29  5 4 3 

30  5 5 4 

31  4 3 3 

32  4 4 3 

33  7 4 3 

34  4 4 5 

33  4 4 i 

34  5 a 3 

37  3 3 3 

33  3 3 3 

39  3 5 5 

40  3 3 3 

41  5 3 5 

42  5 5 3 

43  3 5 5 

44  3 5 5 

43  5 5 3 

44  5 5 5 

47  3 3 3 

48  3 5 3 

49  3 5 5 

30  3 3 5 

31  3 5 5 

32  3 3 3 

33  3 3 5 

34  5 5 5 

35  3 3 5 

54  3 5 5 

57  3 5 5 

58  5 3 5 

59  3 3 3 

40  3 3 5 

41  5 5 3 

42  5 5 5 

43  3 5 5 

44  5 5 5 


•3000E»05  T3T  VOLTAGE 
REAL  1 NAG  SCALE 
7 3 91011121314151617 

55555555555 


5 3 3 3 5 
5 3 5 3 3 
3 3 3 3 3 

3 5 5 5 5 


3 5 3 5 
3 3 5 5 


5 5 3 5 5 
5 5 5 5 5 
3 5 3 5 3 
5 3 5 5 3 


5 5 5 5 3 
5 3 3 3 3 


5 3 5 5 3 
5 5 5 5 3 
3 3 3 4 3 

5 3 3 4 6 
3 5 3 4 4 

6 3 3 3 7 
5 5 4 3 5 
3 5 4 6 5 

3 3 4 4 j 


5 5 5 5 5 
5 3 3 5 3 
5 3 3 3 3 
3 3 3 3 3 
3 3 5 5 5 
5 3 5 5 3 
3 3 5 5 3 
3 3 3 3 3 
5 5 5 5 5 
3 3 3 5 3 
5 3 5 3 3 
3 3 5 3 3 
5 5 3 3 3 
3 5 5 5 3 
5 3 5 3 3 
5 3 3 5 5 


5 3 5 5 3 
5 3 3 5 3 
3 3 5 3 3 
3 5 3 5 3 
5 3 5 3 3 


FACTOR’ 
ISl 9202T ; 
5 5 5 5 
5 5 5 5 
3 5 3 5 


83E-01  0. 

. I OE‘0.1 

222324  32  2 


THIS  PAGE  IS  BEST  QUALlTV  PRACTICABLE 
ERl'M  COri  EUKNISHiD  TO  HOC 


5 5 5 3 
3 5 3 3 
5 5 5 3 
5 5 3 5 
5 3 5 3 
5 5 3 3 
3 5 5 5 
5 5 3 5 
5 5 5 3 
3 3 5 3 
3 5 5 5 
5 3 5 3 


5 3 5 5 
5 3 5 3 
5 3 3 5 
5 5 5 5 


5 5 5 5 
3 5 5 5 
3 5 5 5 
3 3 5 5 
5 5 5 3 
3 3 5 5 
3 5 5 5 
3 3 5 3 

3 3 3 3 

5 3 5 6 

6 5 5 4 

4 6 5 6 
4 3 5 5 
4 6 6 3 

3 3 3 6 

4 5 5 4 


6 3 3 5 
4»65 


5 5 3 5 
3 3 3 5 
5 3 5 5 
5 5 3 3 
5 5 3 3 
3 5 5 3 
5 5 3 5 
5 5 5 3 
3 5 5 5 


3 5 3 6 5 

5 5 6 4 6 
5 5 5 6 3 
5 5 5 5 4 


3 5 5 3 
3 5 3 5 
5 5 5 5 
5 3 3 3 
3 5 3 3 
3 5 3 3 
3 3 3 3 
5 3 5 3 
5 3 3 3 
5 5 5 3 
5 5 5 3 


3 5 3 3 3 
3 3 5 5 5 
5 5 5 5 3 
3 3 5 5 3 
3 3 3 5 3 
5 3 3 5 3 
3 3 5 5 3 
3 5 5 3 3 
5 5 5 5 3 
3 3 5 3 3 
3 3 5 5 3 


INTERVAL  NC.  1 
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Section  5 


CONCLUSIONS 

The  initial  development  of  the  Infrared  Search  and  Track  System 
Evaluation  Model  is  now  complete;  however,  the  refinements  and  modifications 
required  to  make  the  code  more  flexible  and  enable  it  to  reach  its  full  po- 
tential is  far  from  complete.  IRSTEM  contains  a number  of  unique  features, 
which  are  expected  to  make  it  a valuable  tool  for  analyzing  electro-optical 
sensors.  These  features  include  the  detailed,  determinsitic  background 
description  which  obviates  the  need  for  simplifying  assumptions  about  the 
background  sta*’  ;cs.  The  cloud  reflection  and  radiation  computations , 

...e  inclusion  or  _un  c'  fits  and  the  ability  to  input  detailed  sensor  charac- 
teristics are  also  expected  to  be  useful  analysis  features.  . tion,  tne 

separation  of  spectrally-  and  spatially-dependent  effects,  is  expected  to  fa- 
cilitate future  IRSTEM  modifications.  The  overall  goal  of  this  code  is  to 
produce  an  accurate  time-history  of  the  voltage  response  from  each  detector 
so  that  analysts  can  evaluate  different  signal  procession  teonniq  Sucn 
analyses  should  ‘d;et  the  a”  ‘rational  performance  o*  - .r.c  1 • rs 

better  tnan  can  now  oe  done  vitn  tne  statistical  oackgrouiid  tnresn-.j- 
crossing  processing  analysis. 

The  code  described  in  this  report  is  now  operating  but  has  undergone 
only  limited  testing  to  date.  It  is  expected  that  a series  of  modifications 
and  corrections  will  be  required  to  fully  debug  the  program  and  bring  it  to 
its  full  potential.  In  addition,  there  are  a number  of  limitations  associa- 
ted with  the  present  version  of  IRSTEM  which  should  be  corrected  as  time 
allows.  One  of  the  major  limitations  is  the  small  field-of-view  which  is 
necessitated  by  the  tradeoff  between  the  reguired  sensor  resolution  and 
computer  core  limitations.  It  is  felt  that  the  basic  64  by  64  array  of 
0.2  mr  pixels  can  provide  statistically  meaningful  results  without  over- 
loading the  computer.  If  one  wishes  to  increase  the  field-of-view,  this 
can  be  done  by  inputting  a larger  array  (IMAX,  JMAX),  incorporating  larger 
pixels  or  imposing  a moving  background.  The  moving  background  approach 
could  follow  the  target  and  access  information  from  a data  tape  or  simply, 
fold  the  background  back  upon  itself  as  the  target  moves.  The  moving  back- 
ground would  allow  the  user  to  have  an  effective  larger  background  without 
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increasing  the  core  size  or  sacrificing  sensor  resolution. 

Even  with  the  size  background  now  modelled  in  IRSTEM.  the  code 
may  be  too  large  for  some  users.  Work  is  underway  to  reduce  the  code 

s _e  by  shortening  the  LCWEV  subroutine  and  reducing  the  sun  glint  calcula- 
tions. Anotner  option  is  to  split  the  code  into  two  or  mere  pieces.  Referr- 
ing back  to  Figure  23,  one  can  see  that  there  is  a natural  break  Doint 
between  the  one-time-only  calculations  and  the  major  time  loop.  In  addition, 
if  the  target  signature  and  glint  calculations  are  moved  to  the  one-time 
only  calculations,  the  calls  to  LOWEM  can  be  eliminated  in  the  time-loop, 
and  portions  of  the  program  can  be  overlaid  to  reduce  the  core  size 
requirements. 

Another  limitation  is  the  fact  that  the  code  models  only  radially 
inbound  missiles.  The  graybodv  signature  calcu'ation  used  for  these  targets 
is  a simplistic  one,  and  the  computations  terminate  nen  the  target  exceeds 
the  pixel  size.  Finally,  jitter  and  turbulence  have  not  been  incorporated 
into  the  code  as  of  this  time. 

Work  is  underway  to  expand  the  capabilities  of  IRSTEM  and  eliminate 
many  of  the  limitations  identified  above.  However,  even  when  these  modifica- 
tions are  completed,  there  is  still  a large  number  - :r:;cns  available 
for  increasing  the  flexibility  and  utility  of  IRSTEM.  Fm-  example,  one 
may  wish  to  incorporate  generalized  trajectory  and  signature  computations 
which  are  time-varying  and  allow  the  user  to  model  target  maneuvers  and 
missile  burnout.  Another  growth  option  to  pred’ :*■  alarm  rates. 

The  present  code  is  designed  to  analyze  the  probability  of  tarqet  detec- 
tion against  a variety  of  backgrounds  but  it  is  not  optimum  for  analyzing 
false  alarms  because  of  its  deterministic  approach.  There  are  techniques, 
however,  which  can  be  used  to  obtain  statistical  information  on  the  false 
alarms  using  the  random  noise  portion  of  the  program.  Finally,  even 
though  the  present  code  is  structured  to  model  the  wide  field-of-view 
surveillance  systems,  the  basic  methodology  is  also  applicable  to  FURS, 
reticle  seekers  and  non-scanning  sensors. 
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At  present,  it  appears  that  IRST"M  testing  should  continue  using 
existing  hardware  and  field  test  data  whenever  available.  In  addition, 
it  is  now  ready  to  predict  the  performance  of  advanced  concepts.  Performing 
such  analyses  will  not  only  uncover  subtle  coding  errors  but  also  will  lead 
to  logic  modifications  which  will  improve  tne  accuracy  and  utility  of 
IRSTEM. 
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